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NUCLEAR  REACTOR  SIMULATOR 


w 


Preface 


The  purpose  of  this  report  is  to  describe  the  work  done  at  Oak  Ridge 
National  Laboratory  in  the  development  of  a  nuclear  reactor  simulator. 

This  work  is  part  of  the  total  program  of  the  Power  Pile  Division  and  was 
carried  on  by  the  Control  Section.  No  one  person  in  this  section  spent  full 
time  in  this  development  work.  Limited  personnel  made  it  necessary  to 
stop  this  work  from  time  to  time  to  design  electrical  control  circuits  for 
the  particular  reactors  being  studied  by  the  division.  The  advantages  to  be 
gained  from  having  available  a  simulator  were  sufficient  to  warrant  direct¬ 
ing  the  persons  involved  to  put  as  much  effort  as  possible  into  its  develop¬ 
ment. 


A  great  deal  has  been  learned  about  the  possibilities  of  such  a 
simulator  and  the  difficulties  which  are  encountered  in  its  design.  In  this 
report  no  attempt  is  made  to  draw  any  definite  conclusions.  It  is  felt  that, 
because  of  changes  in  personnel  and  the  fact  that  the  reactor  development 
program  has  been  relocated,  all  the  significant  information  pertaining  to 
the  simulator  should  be  gathered  together. 

This  report  includes  the  following: 

Introduction,  Page  7  ,  is  a  general  picture  of  the  ideas  behind  the 
simulator  and  its  uses. 

Section  A,  Page  14,  entitled  Network  Development,  covers  the  der¬ 
ivation  of  formulae  for  the  relations  between  reactor  constants  and  elec¬ 
trical  network  constants. 

Section  B,  Page  28  ,  entitled  Construction  of  Experimental  Network, 
describes  the  networks  used  to  determine  the  feasibility  of  the  simulator. 
Results  of  test  runs  are  given  and  compared  to  calculated  values. 

Section  C,  Page  45,  entitled  Current  Feeding  Equipment,  gives  a 
description  of  the  methods  used  to  feed  electrical  current,  which  simulates 
the  birth  of  neutrons  in  the  reactor,  into  the  electrical  network. 

Section  D,  Page  65,  is  entitled  Experiment  and  Development  Work, 
As  the  tests  described  above  were  conducted  several  ideas  developed  which 
offered  solutions  to  the  various  difficulties  encountered.  Some  of  these 
were  tested  and  are  described  in  this  section. 


Appendix  I,  Page  90,  is  entitled  Transmission  Line  Analogy,  In 
this  section  an  attempt  is  made  to  present  corrective  formulae  which  would 
make  the  results  of  measurements  on  a  lumped  network  exactly  equal  to 
the  equations  for  a  distributed  network  and  therefore  equivalent  to  a 
homogeneous  reactor. 
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SYMBOLS  AND  DEFINITIONS  OF  TERMS 


General  symbols  used  for  reactor  constants  and  quantities  through¬ 
out  this  report  are: 

k  =  number  of  thermal  energy  neutrons  produced  per  thermal 
neutron  absorbed  in  an  infinite  core. 

Ztth  “  transport  cross  section  for  thermal  energy  neutrons,  (cm-1). 

=  transport  cross  section  for  high  energy  neutrons,  (cm*1), 

2ath  =  absorption  cross  section  for  thermal  energy  neutrons,  (cm"1) 

2af  =  absorption  cross  section  for  high  energy  neutrons,  (cm"1). 

$£  =  high  energy  neutrons  per  sq  cm  per  second. 

$  =  thermal  energy  neutrons  per  sq  cm  per  second. 


32tth2ath  6 


of  the  mean  square  distance  a  neutron  travels 


from  point  at  which  it  becomes  thermal  to  point  at  which  it 
is  absorbed. 

T  =  Yv - v —  =  "Fermi  age”  of  neutrons,  of  the  mean  square 

32jtf  af  6 

distance  a  neutron  travels  between  point  of  fission  and  point 
at  which  it  becomes  a  thermal  energy  neutron. 

The  terms  used  to  describe  the  electrical  system  are  those  generally 
accepted  in  electrical  engineering.  However,  the  terms  are  somewhat  more 
specific  and  are  set  out  here  for  clarity. 

Network  is  the  term  applied  to  that  portion  of  the  simulator  consisting 
of  resistances  interconnected  in  a  mesh  configuration. 

Simulator  is  the  term  applied  to  the  entire  system  which  includes 
the  network  and  the  current  feeding  equipment. 

Node  Point  is  the  term  describing  those  points  in  the  resistor  net¬ 
work  where  resistors  are  joined  together. 


Lattice  Unit  is  the  term  used  to  describe  that  portion  of  the  network 
of  resistances  which  contain  a  node  point.  In  particular  refer  to  fig,  A-l. 


The  lattice  unit  here  consists  of  the  four  resistances  r/Z  and  the  resistor 
Rg.  In  the  network  this  point  would  be  surrounded  by  four  similar  lattice 
units  a  distance  h  away  from  x,  y,  which  would  contain  the  other  half  of  R 
shown  here.  Each  of  these  four  lattice  units  also  contain  the  value  Rg. 
Fig.  A-Z  is  a  typical  lattice  unit  for  a  three  dimensional  network.  In  this 
figure  two  lattice  units  similar  to  fig.  A-l  have  been  placed  on  top  of  one 
another  and  a  lattice  unit  refers  to  this  combination.  Again,  each  lattice 
unit  contains  the  value  R  and  Rg  and  six  of  each  of  the  values  Rg/2  and 

Feed  Point  is  the  term  used  to  indicate  the  node  points  at  which 
current  is  fed  into  the  network.  Feed  points  occur  only  in  the  network 
which  represents  reactor  region. 
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NUCLEAR  REACTOR  SIMULATOR 
Introduction 


The  simulator  is  an  electrical'analogue  computer.  The  idea  of 
solving  differential  equations  by  the  use  of  electrical  models  has  become 
a  widely  accepted  technique.  If  due  consideration  has  been  given  to  keeping 
the  model  similar  to  the  physical  phenomena  described  by  the  differential 
equations,  the  practical  engineer  has  at  his  disposal  a  means  of  studying 
the  physical  phenomena  without  continual  reference  to  the  differential 
equations  involved.  The  simulator  developed  at  ORNL,  through  the  combined 
efforts  of  E.  J.  Wade  and  J.  W.  Simpson  of  the  Power  Pile  Division,  meets 
this  requirement.  It  consists  of  a  network  of  electrical  conducting  ele¬ 
ments  and  suitable  sources  of  current  and  current  sinks.  The  design  of 
this  network  is  such  that  when  the  proper  electrical  resistance  values  are 
assembled  and  the  current  sources  and  sinks  are  properly  adjusted,  the 
steady  state  voltage  distribution  on  the  network  is  the  same  as  the  steady 
state  flux  distribution  in  the  reactor  and  the  current  sources  are  a  measure 
of  the  critical  mass  of  the  reactor. 

In  elementary  diffusion  theory,  the  neutrons  are  assumed  to  obey 
the  diffusion  equation  (a)  for  any  volume. 

3  dt 

where 

A 


3> 


Za  =  absorption  constant  for  neutrons. 

Q  =  neutrons  produced  per  cu  cm  per  sec  which  may  be  a 
function  of  their  position  in  the  reactor. 

n  =  number  of  neutrons  per  cu  cm. 

t  =  time. 

V  =  Laplacian  =  divergence  of  the  gradient. 

(a)  AECD-2201 

6:76 


=  ^a  = 

At 

3  " 


number  of  neutrons  absorbed  per  cu  cm  per  sec. 
diffusion  constant  for  the  reactor  material. 


=  (nv)  =  neutron  flux  -  neutrons/cu  cm  x  cm  per  sec  =  neutrons 
per  sq  cm. 


This  equation  shows  that  the  rate  of  change  of  neutrons  for  a  par¬ 
ticular  volume  of  core  material  equals  the  net  neutrons  that  diffuse  into 
that  volume  less  those  that  are  absorbed  there  plus  those  produced  in  the 
same  volume. 


The  net  neutron  current  density  is  given  by  equation  (b), 


The  diffusion  of  electrical  current  in  a  small  volume  of  electrical  conductor 
obeys  the  diffusion  equation  (c), 


+  I 


dQ 

dt 


in  which 

V  =  electrical  potential. 

I  =  electrical  current  entering  the  conductor. 

Q  =  electrical  charge  per  unit  volume, 
t  =  time. 

V2  =  Laplacian. 

R  =  resistance  of  the  conductor. 

C  =  the  sink  or  loss  of  current  per  sec  for  that  volume  of  conductor. 
The  current  flowing  in  the  conductor  is  given  by  equation  (d) 

V  V 

I  =  — R- 

The  similarity  of  equations  (a)  and  (c)  as  well  as  (b)  and  (d)  can  readily  be 
seen. 

V  is  analogous  to  (nv) 
i  is  analogous  to  J 

Q  is  analogous  to  n 

A  volume  of  electrical  conductor  in  which  there  is  a  sink  C  and  a 
source  I,  will  serve  as  an  electrical  analogue  for  the  differential  equation 
of  neutron  diffusion.  Consideration  is  restricted  to  only  the  steady  state 

(b)  AECD-2201 
6:76 
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flux  condition  and  therefore  in  equation  (a)  and  (c)  the  right  hand  side  is 
zero. 


* 


w 


3 

VZV 

R 


V2<I>  -  3>2a 
-  C  +  I  =  o 


+ 


Q  = 


O 


(e) 


(f) 


Practical  aspects  of  building  a  conductor  such  as  have  been  de¬ 
scribed  must  now  be  considered.  It  is  not  difficult  to  build  a  conductor  of 
some  material  and  introduce  a  source  or  several  sources  of  current  into 
it  through  small  insulated  conductors.  However,  it  is  quite  impossible  to 
construct  an  absolutely  homogeneous  distribution  of  sources  in  the  conductor. 
Therefore,  a  finite  spacing  of  sources  must  suffice.  Since  the  sources  must 
be  spaced  at  finite  intervals  the  finite  volume  of  conductor,  which  exists 
between  each  source  may  as  well  be  replaced  with  a  conductor  which  has 
the  same  resistance  as  that  volume  of  conductor.  If,  for  simplicity,  the 
sources  are  spaced  evenly  throughout  the  volume,  the  same  resistor  be¬ 
tween  each  source  would  be  used.  So  far  it  has  been  implied  that  the  cur¬ 
rent  is  being  removed  from  all  points  on  the  surface  of  the  volume  of  con¬ 
ductor  and  that  these  points  are  not  necessarily  at  the  same  potential.  The 
same  is  true  for  the  network  of  lumped  conductors.  The  current  is  removed 
from  all  resistors  which  reach  to  the  boundary  of  the  volume.  For  the 
balance  of  this  discussion  the  distance  between  the  two  sources  is  considered 
to  be  a  unit  distance  (assuming  evenly  spaced  sources)  and  a  volume  of  all 
unit  dimensions  is  called  a  unit  volume. 


For  convenience  attention  is  confined  to  a  volume  of  conductor  of 
unit  thickness  which  contains  all  the  sources  in  one  flat  plane  of  the  total 
conductor  volume.  Here  it  is  implied  that  the  components  of  current  flow 
out  of  this  unit  thickness  volume,  perpendicular  to  the  plane  of  the  sources, 
are  known.  A  network  to  represent  such  a  volume  is  shown  in  fig.  1  where 
equal  spacing  of  the  sources  is  again  used  and  therefore  all  the  conductors 
are  of  equal  resistance  R  ohms.  At  each  of  the  node  points  (where  resistors 
join)  the  current  source  I  is  .introduced.  This  current  flows  through  the 
network  eventually  reaching  ground  to  complete  the  electrical  circuit.  This 
current  I  can  be  any  quantity  and  all  the  feeds  are  not  necessarily  equal. 

For  the  analogy  this  current  should  be  proportional  to  the  voltage  which 
exists  at  each  node  point  because,  in  the  reactor,  the  number  of  neutrons 
produced  by  fission  is  proportional  to  the  thermal  flux  present  at  any  point 
in  the  reactor,  where  fissions  occur.  To  complete  the  analogy  a  sink  term 
C  is  needed,  which  for  a  perfect  analogy,  requires  that  the  electrical  current 
be  destroyed  at  all  points  in  the  volume  of  conducting  material  in  an  amount 
which  is  everywhere  proportional  to  the  neutron  flux  present.  Here  again 
practical  considerations  require  that  the  sink  C,  for  the  units  of  volume  be 
lumped  at  the  node  points  of  the  network.  Since  this  sink  is  proportional  to 
the  voltage,  at  each  node  point,  it  is  possible  to  take  advantage  of  the  linear 


P,P,P.A-l4i</  Fig. 


relation  of  Ohm’s  law 

I  =  e/r 

and  place  a  conductor,  from  each  node  point  to  ground,  which  has  a  resist¬ 
ance  inversely  proportional  to  the  absorption  constant  for  neutrons.  The 
current  source  requirement  can  be  met  by  a  voltage  controlled  current  sup¬ 
ply  such  as  a  vacuum  tube  operating  in  the  linear  portion  of  its  characteristic 
curve.  With  a  suitably  controlled  source  in  place,  the  analogy  for  elementary 
diffusion  theory  is  complete. 

In  a  thermal -neutron  nuclear  reactor,  neutrons  exist  at  all  energies 
between  fission  energy  and  thermal  energy,  and  the  '‘constants39  of  equation 
(e)  are  functions  of  neutron  energy.  The  atwo-group  theoryS3*is  used  exten¬ 
sively  in  practical  reactor  calculations  as  an  approximation  to  this  situation. 
Here  it  is  assumed  that  the  neutrons  exist  in  two  energy  groups,  fast  or 
high  energy  and  slow  or  thermal  energy.  For  a  steady  state  flux  of  either 
energy  level,  an  equation  like  e  applies.  In  the  case  of  fast  neutrons  the  A 
term  represents  slowing  down  while  for  thermal  neutrons  it  continues  to  be 
absorption.  The  electrical  analogy  requires  two  similar  networks  of  the 
type  described.  Each  must  have  the  same  lattice  construction,  that  is, 
equally  spaced  node  points.  Thermal  neutrons  are  produced  by  the  slowing 
down  of  fast  neutrons.  The  current  of  the  sink  term  of  the  network  repre¬ 
senting  fast  neutrons  must  be  made  to  flow  into  the  corresponding  node 
points  of  the  thermal  network.  This  suggests  that  the  conductor  which 
carries  the  sink  current  away  from  the  fast  network  should  connect  directly 
between  the  corresponding  node  points  of  the  fast  and  thermal  networks 
instead  of  to  ground.  In  two-group  theory  it  is  assumed  that  all  fissions 
occur  because  a  thermal  neutron  flux  is  present  at  any  point  at  which  fis¬ 
sion  can  occur.  For  this  analogy,  then,  the  current  source  I,  for  the  fast 
network  is  controlled  with  the  voltage  of  the  thermal  network.  Thus,  a 
single  network  can  be  constructed  to  simulate  all  the  elements  of  two-group 
diffusion  theory. 

A  typical  lattice  unit  for  this  network  is  shown  in  fig,  A-l  of  Section 
A.  In  this  figure  the  fast  neutron  diffusion  network  resistors  are  denoted 
by  the  subscripts  0  and  the  thermal  network  by  6  .  R  ^  is  the  slowing  down 
resistor  and  Rg  is  the  absorption  resistor.  The  source  current  I  is  repre¬ 
sented  by  a  box  which  is  an  electronic  amplifier  whose  output  current  I  is 
equal  to  A  x  Vg  .  To  operate  the  simulator  a  particular  value  of  A  must  be 
found  which  will  establish  constant  voltages  on  all  points  in  the  network. 

If  the  voltage  V Q  were  zero  the  source  I  would  be  zero.  In  order  to  start 
the  cycle  of  production  of  current,  a  sourve  Y Q  must  be  applied  to  the  net¬ 
work  just  as  in  the  reactor  a  source  of  neutrons  must  be  present  in  order 
to  start  the  chain  reaction.  Once  the  cycle  is  started  the  source  can  be 
removed  and  the  cycle  simulating  the  chain  reaction  will  continue. 


*Soodak,  H.,  AECD  Publication  2201 


An  ideal  analogue  would  simulate  the  continuous  slowing  down  proc¬ 
ess  which  occurs  for  the  neutrons  in  the  chain  reactor.  Actually  what  is 
wanted  is  the  simulation  of  many  groups  of  neutron  energies.  From  the 
discussion  of  the  two-group  simulator  it  is  evident  that  a  voltage  difference 
must  exist  such  that  the  0  (fast)  network  voltage  is  higher  than  the  9 
(thermal)  network  voltage.  The  larger  this  voltage  difference,  the  more 
nearly  the  slowing  down  of  neutrons  can  be  simulated  in  the  sink  conductor 
of  the  fast  network.  In  the  simulator  of  this  report  (Section  B)  the  0  net¬ 
work  peak  voltage  was  set  at  150  volts  and  the  9  network  peak  at  3  volts, 
a  ratio  of  50:1.  Obviously  an  attempt  to  simulate  three  groups  with  the 
same  ratios  results  in  ratios  of  2500:50-1,  a  peak  of  7500  volts  for  the 
highest  energy  for  a  peak  of  3  volts  on  the  lowest  energy  network.  This 
is  not  practical  but,  fortunately,  the  50:1  ratio  between  voltages  is  not 
necessary  and  by  using  a  smaller  ratio  of,  say  10:1,  at  least  three  energy 
groups  could  be  simulated.  At  ORNL  the  work  was  confined  to  two  energy 
groups . 


In  Section  A  of  this  report,  the  relations  between  the  network  param¬ 
eters  and  reactor  constants  are  developed.  In  Section  B  the  construction 
of  two  simulators,  each  for  a  different  reactor,  is  described. 

A  chain  reactor  consists  of  a  core  and  a  reflector  region.  So  far 
only  a  region  in  which  sources  of  neutrons  are  present,  which  is  the  case 
in  the  core,  has  been  considered.  The  slowing  down,  absorption  and  dif¬ 
fusion  which  occur  in  the  reflector  region  can  be  simulated  by  a  network 
of  the  same  form  as  the  core  with  the  source  current  (i)  zero.  The  source 
of  thermal  neutrons  due  to  slowing  down  is  still  present  and  the  absorption 
of  thermal  neutrons  is  much  less  because  of  the  absence  of  fissionable 
material.  For  all  practical  purposes,  all  that  is  necessary  in  order  to 
simulate  the  reflector  region  is  to  assign  average  values  of  material  con¬ 
stants  to  that  region  and  to  adjust  the  network  resistance  values  to  corre¬ 
spond.  This  is  covered  more  completely  in  the  designs  of  Section  B.  The 
neutron  flux  is  considered  to  be  zero  at  the  outside  boundary  of  the  reflector 
and  therefore  the  network  is  grounded  at  that  point. 

Simulating  the  insertion  of  a  control  rod  into  the  core  is  a  matter 
of  imposing  upon  the  simulator  network  the  known  effect  of  the  rod  upon 
the  neutrons.  If  the  rod  will  absorb  all  the  thermal  neutrons  in  a  particular 
reactor  region  the  thermal  network  can  be  grounded  at  those  points  which 
correspond  to  the  effective  rod  dimensions.  The  practical  considerations 
of  this  operation  are  covered  more  fully  in  Section  B  of  this  report. 

With  a  simulator  of  this  kind  available,  the  determination  of  critical 
mass  and  the  flux  distribution  for  complex  boundary  condition  problems  is 
greatly  simplified.  Analytical  solution  of  the  two-group  theory  differential 
equations  results,  for  nearly  all  boundary  problems,  in  a  tedious  trial  and 
error  process.  The  simulator  offers  a  means  of  performing  this  process 


rapidly  through  the  feed-back  action  of  the  electronic  amplifier.  Examples 
of  these  complex  problems  are:  the  insertion  of  control  rods  along  the 
radius  or  a  chord  of  a  cylindrical  core,  irregular  control  rod  spacings  and 
the  case  where  the  material  constants  vary  as  a  function  of  the  location  of 
the  material  in  the  reactor. 
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SECTION  A 

NETWORK  DEVELOPMENT 


A.  Network  Theory 

Partial  differential  equations  such  as  those  used  in  two-group  theory- 
calculations  can  be  solved  by  the  use  of  suitable  electrical  networks. 

Consider  the  point  x,  y  in  a  network  mesh  as  shown  in  fig.  A-l  (b) 

(This  is  one  point  in  the  network  of  fig.  A-l  (a). 

The  electrical  current  equation  at  x,  y  is: 

V (x+h,y)  -  V (x,y)  ,  V(x-h,y)  -  V(x,y)  ,  V(x,y+h)  -  V(x,y)  .  V(x,y-h)  -  V(x,y)  =  Vfcjrj 

- r - +  r  +  R  R  Rg 


or; 

R 

[V (x+h,y)  +  V (x-h,y)  -  2V(x,y)]  +  [V(x,y+h)  +  V(x,y-h)  -  2V(x,y)]  -  V(x,y)  — 

o 


Using  Taylor's  expansion* 


V  (x+h,y)  =  V(x,y)  +  -yj 


dV(x,y)  hf  d  ZV (x,y) 
3x  2!  ax2 


and 


V  (x-h,y)  =  V  (x,y) 


h  av(x,y)  h2d2V(x,y)  _  _  _ 

i!  ax  2!  dxz 


adding  and  rearranging  terms 

3  2V (x,y)  _  V (x+h,y)  V(x-h,y)  2V(x,y) 

3x2  h2 


V 


*In  this  development  a  lattice  length  h  =  1  could  be  conveniently  used,  and 
the  variables  x,  y  and  z  could  be  measured  in  lattice  units.  In  Section 
B,  however,  different  values  of  length  for  h  were  used  and  it  was  considered 
more  convenient  to  let  h  remain  as  a  parameter  in  the  equation. 
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likewise 


32V(x,y)  _  V(x,y+h)  +  V(x,y-h)  -  2V(x,y) 


The  network  equation  becomes  very  nearly: 

aZV(x,y)  d2V(x,y)  _  V(x,y)  JR 
b  x2  d  y2  h2  Rg 

The  solution  of  this  equation  V(xy)  can  be  read  with  a  voltmeter 
from  the  network  point  (x,  y)  to  ground. 


B.  Network  Theory  Applied  to  Calculations  for  Reactor 

This  fundamental  idea  can  be  extended  into  three  directions  x,  y 
and  z  and  the  resulting  network  equation  will  have  aV2V  term  where  the 
reactor  equation  contains  V 2 The  network  equation  would  be: 

v  2V(x,y,z)  .  rR 

h  Rg 

In  two-group  theory  the  neutrons  are  assumed  to  exist  only  in  two 
energy  levels,  fast  or  high  energy,  and  thermal  or  low  energy.  In  either 
energy  level  the  general  form  of  the  balance  equation  at  steady  state  condi 
tions  of  flux  is: 

-  Leakage  -  absorption  +  production  =  0 


<‘p> 


3  ^^th  ^  ath^  th  +  ^af 


where 


3  f 


number  of  fast  neutrons  diffusing  into  a  cu  cm  per  sec 


LI1_2 

~  number  of  thermal  neutrons  diffusing  into  a  cu  cm  per  sec 


2af  $£  =  number  of  fast  neutrons  becoming  thermal/cu  cm/sec. 
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^ath^th  =  number  of  thermal  neutrons  absorbed/ cu  cm/sec. 

$  =  number  of  fast  neutrons  produced/cu  cm/sec. 

ath  th 

In  the  electrical  network  a  current  1  amp  can  then  be  introduced  at 
the  node  point  and  the  balance  equation  for  current  will  become: 

*^V2V(xyz)  -  +  I  -0 


or; 


H72V(xyz)-V(xyltlz0 
R  hRg  h 

which  is  of  the  same  form  as  (Ip)  or  (2p)  above. 


For  each  neutron  energy  group  an  electrical  resistor  network  having 
the  above  equation  can  be  constructed.  For  a  simultaneous  solution  of  these 
equations,  the  current  I  at  each  node  point  in  each  network  must  be  made 
proportional  to  the  voltage  V  of  the  corresponding  node  point  in  the  other 
network. 

* 

Rather  than  construct  separate  networks  to  meet  these  requirements, 
a  single  network  with  suitable  resistance  elements  to  represent  all  of  the  V 

terms  in  both  energy  group  equations  can  be  arranged.  The  network  equations 
can  then  be  written  and  the  coefficients  of  these  equations  equated  to  the 
coefficients  of  the  reactor  equations. 

A  network  of  resistances  to  meet  these  requirements  is  shown  in 
fig.  A-2. 


C.  Calculation  of  Network  Constants 


The  current  equation  at  V,  \  is: 

\x,y  ,z) 

(Equation  In) 

V<j>(x+a,y,z)  ~  Vrf,(x,y,z)  V^(x-a,y,z)  -  V(ft(x,y,z)  + 

0x  </>x 

V^x.y+b.z)  -  V<j(x,y,z)  ^  V^(x/y-b,z)  -  V^(x,y,z) 

—  - R  R  , 

0y  0y 

V<Mx>y,z+c)  -  v^(x,y,z)  V,* (x,y,z~c)  -  V^(x;y,z)  +  V^Vg  +  AVq  =  0 


R 


<p  z 


R 


R 


7 
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or  collecting  terms 
V0(x+a,y,z)  +  V0(x-a,y,z)  -  2V0  (x,y,z) 


V0(x,y+b,z)  +  V0(x,y-b,z)  -  2V0(x,y,z) 


V0(x,y,z+c)  +  V0(x,y,z-c)  -  2V0(x,y,z)  V0 


+  +  V0  [  A  -  -g-  ]  =  0 


Reactor  equation  for  fast  flux  (nv)f  =  Of 


Xtf  a2 


&  "2  af$f  +  kSath^th  =  0 


By  Taylor  Expansion  (using  a,  b,  and  c  to  correspond  to  lattice  unit  length 
in  x,  y,  and  z  directions) 

Of(x+a,y,z)  =  Of(x,y,z)  +-~j  yj^(x,y,z)  +  J\  g  (x,y,z)  +  .  .  .  . 

Of(x-a,y,z)  =  €>f(x,y,z)  (x,y,z)  +  ”  '  '  '  ’ 


Adding 

d2$f  (x,y,z)  Of(x+a,y,z)  +  Of(x-a,y,z)  -  20f(x,y,z) 

dP  P 

Similarly,  terms  for 


3  2  Of(x,y  ,z) 


32$f  (x,y,z) 

—  a? — 


can  be  obtained. 


Substituting  these  expressions  in  lp  we  get:  (after  multiplying 
through  by  a,  b  and  c) 

Equation  Ip 


tf  be 
3  a 


[  Of(x+a,y ,z)  +  Of(x-a,y,z)  -  2<I> f(x,y ,z)] 


tf.  ac 

T’T 


[ <l>f  (x,y+b,z)  +  Of  (x,y-b,z)  -  20f(x,y,z)]  + 
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Xtf  ab 


3  c 


[<E>  (x,y,z+c)  +  4>  (x,y,z-c)  -  20  (x,y,z)]  - 


abc  2  $£(x,y,z)  +  abck  2a£h  <|>th(x,y,z)  =  0 

Equation  In  and  Ip'  will  be  equal  if: 

•f  =  V<t> 


(4-1) 

(4-2) 

(4-3) 

(5) 

(6) 


0 

th 

R 

II 

a 

0x 

*tf 

#  be 

RJ 

3 

b 

<t>  y 

^-tf 

ac 

R 

3 

c 

0  Z 

Xtf 

’  ""ab 

R7  2  „  abc 


—  +  A  =  k  2  abc 
R-y  ath 


Likewise,  the  equation  for  the  9  network  is  as  follows:  (Using  a,  b  and  c 
as  before). 

(Equation  2n) 

Vq (x+a,y,z)  +  V^x-a.y.z)  -  2Vfi(x,y,z) 

Rex 

Vfi(x,ytb,z)  +  V o(x,y-b,z)  -  ZVe(x,y,z) 


Vfi(x,y,z+c)  +  V  fl(x,y,z-c)  -  2Ve(x/y,z)  _  ,  (  f  JL  x  J_1  x  =  n 

The  reaction  equation  for  thermal  flux  ($)£h  infinite  difference  form  is: 

(Equation  2p') 

— ^  [®th(x+a,y,z)  +  <E>th(x-a,y,z)  -  20th(x,y,z)]  + 


^•tth  ac 


.  —  [4>  (x,y+t>,z)  +  <1>xi.(x*y"b»z)  "z$+x(x»y»z))  + 


3  b  1  th 
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— y1  •  [  't’th(x.y.z+c)  +  ^(x.y.z-c)  -  2Hth(x,y,z)] 

abc  Sath  '*,th<x’y'z)  +  abc  Zaf  Vx'y,!!)  =  ° 

Here  again  Zn  will  be  similar  to  Zp'  if: 


V6^th 

R  -  3  a 

^  "  \th  ‘  ^ 

R  -_L  _L 
0y  xtth’  ac 

R  -  3  C 

rr  ■  “sh" 

tth 

1  1  rr 

Rg  Ry  7ath 


(7) 

(8-1) 

(8-2) 

(8-3) 

(9) 


Having  established  the  above  relations  the  ratios  between  network 
resistance  values  can  be  set  up. 


D.  Network  Resistor  Ratios 


R  .  3a*  ^  di 

- 2af  =  — 

R7  Xtf  T 

K.  3b2  _  b2 

*tf  af  T 

^L=2l‘z  f-i! 

R7  Xt£  af  T 

From  4  -  1,  Z,  3  and  8  -  1,  Z  and  3 
R  j  R .  R .  A. 

0x  0y  0z  _  t  ,th 


@x  ^0v  ^0z  A.  tf 


From  8-1  and  7 


r  i  ii 

to 

CD 

kgV* 

R0X 

1  — 

Rg 

*7  ~L2 

'ath  L,2 


(14) 
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From  13 


CD 

X 

_  «  Xtf 

K0X  • 

X  tth 

R  ~ 

2 

-  JL- 

Rcfex 

0x 

"  L2 

R7 

R0x 

_  a2 

"l  - 

Rg 

‘  L2 

T 

1 


T  Xtth 


(14-1) 


Likewise: 


"0  y  _  b 


L2  Xtf 
T  X  ,+1 


(14-2) 


^ r, 

Rg  L‘  l 


_l2  Xtf 
T  X  tth. 


From  6 

A  =  k 

Substitute  9 


ath  Ry 


A-k[V<K 


k_  k  -  1 
Rg  R7 


The  current  I  was  defined  as  AV. 


A  =  -  (mhos) 

V0 


(14-3) 


In  order  to  supply  this  current  a  current  generator  is  required  which 
can  be  controlled.  Its  output  current  must  always  equal  A Vq  .  An  electronic 
amplifier  will  satisfy  this  requirement.  The  ideal  input-output  relation  for 
this  amplifier  is  a  straight  line,  the  slope  of  which  is  A.  The  slope  A  must 
be  adjustable. 
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E.  Kappa  Factor 

The  network  under  discussion  will  simulate  a  reactor  described  by 
the  two-group  theory  equation  but  the  current  flows  iy,  ig  and  I  of  fig.  A-2 
will  not  represent  exactly  any  reactor  quantities  because  the  0  network  has 
been  placed  on  top  of  the  0  network.  Consider  again  the  basic  network  re¬ 
quirement  that  the  absorption  resistor  (R/y  or  R$)  should  go  to  ground.  It 
is  easy  to  see  that  this  ideal  network  will  be  approximated  by  the  combined 
network  if  the  voltage  V0  of  this  combined  network  is  very  large  compared 
to  V$  ,  for  example,  in  the  order  of  a  50:1  ratio. 

This  means  that  there  must  be  a  ratio  of  about  50:1  between  the 
values  of  Ky  and  Rg  and  between  the  values  of  R0  and  Rg. 

Formula  (l3)(with  x  =  y  =  z)  gives  the  relation 

Atth 
Ttf 

A  factor  k  is  now  introduced  so  that 

R<£  _  k  Xtth 
R0  xtf 

and  this  is  called  a  “kappa  factor.”  This  gives  the  desired  arbitrary  re¬ 
duction  of  the  network  voltage  and  if  ^tth  =  ^tf 


(13a) 


R<f) 

R-e 


R$  =  kR@ 

This  factor  also  affects  formulas  (14)  which  (with  a  =  b  =  c)  becomes 

Re  = a2  r  _  L2  Xtf_1 

Rg  L2  L  T  XtthJ 

and  before  the  “kappa  factor”  was  introduced  this  could  have  been  written 
(by  using  13) 

M  .Rfll 

Rg  L2[_  T  R<*>J 

Now  inserting  the  ratio  of  13a  with  the  “kappa  factor”  introduced  the  result 
is: 

Rg  L2  L  T  ^tth 


(14a) 
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In  formula  15  the  value  of  Rg  which  was  determined  using  k  is  em¬ 
ployed.  Here,  since  Rg  has  been  reduced,  A  for  a  given  k  is  larger.  This 
means  simply  that  any  method  of  introducing  the  current  I  =  A Vq  must 
supply  a  larger  current  per  volt  than  would  have  been  required  without  the 
kappa  factor.” 

When  the  relations  of  formulas  3  through  15  are  built  into  the  elec¬ 
trical  network  the  following  direct  relations  between  reactors  and  network 
quantities  exist. 


Network  Quantity  Reactor  Quantity  Represented 


i^x.  i<£y»  i0z; 

Diffusion  of  fast  neutrons/sec/a.b.c. 

(cm)3  in  the  directions  x,  y  or  z 

iy 

Fast  neutrons/sec/a.b.c.  (cm)3  which 
slow  down  or  become  thermal 

i©x>  *6y>  * 9z 

Diffusion  of  thermal  neutrons/sec/a.b.c. 
(cm)3  in  the  direction  x,  y  or  z 

Thermal  neutrons/sec/a.b.c. (cm)3  which 
cause  fissions 

I  (The  current  entering  the 

Fast  neutrons/sec/a.b.c.  (cm)3  which 

network  node  point) 

result  from  fissions 

A  solution  to  the  network  equations  results  when  V<£  and  VQ  are 
constant  and  for  this  condition  A  for  the  amplifier  must  be  adjusted  to  a 
particular  value.  The  operational  procedure  is  to  start  current  flowing  in 
the  network  by  some  external  means,  and  vary  A  until  V<£  and  VQ  are 
constant.  A  is  then  a  measure  of  k  for  the  reactor  by  equation  15.  The  node 
point  voltages  V 0  and  Vq  give  the  flux  distribution  in  the  reactor. 


F.  Boundary  Conditions 

So  far  only  the  core  has  been  considered  and  no  specific  attention 
has  been  given  to  the  reflector  region.  The  core  is  surrounded  by  a  reflec¬ 
tor  region  and  it  is  assumed  that  no  neutrons  leave  the  reflector.  Three 
boundary  conditions  regularly  used  for  solution  of  the  reactor  equations  must 
be  met  by  the  simulator;  namely,  that  in  the  reflector  the  production  of 
neutrons  is  zero,  that  the  flux  is  continuous  across  the  interface  between  the 
reflector  and  the  reactor  and  that  the  thermal  and  fast  flux  are  both  zero  at 
the  outside  boundary  of  the  reflector.  The  first  condition  is  met  by  not 
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feeding  current  into  the  reflector  node  point.  To  meet  the  third  of  these 
conditions  the  perimeter  of  the  network  which  represents  the  reflector  is 
grounded.  For  the  second  condition  reflector  section  network  must  be 
built  to  contain  the  constants  T»  L2,  and  \th.  for  the  reflector  using  the 
formulas  developed  for  R0,  Rq,  Ry,  and  Rg  and  the  network  for  the  reactor 
is  joined  to  this  reflector  network.  Due  consideration  must  be  given  to 
the  location  of  the  interface  with  respect  to  the  network  node  points. 


If  the  interface  occurs  midway  between  node  points  the  value  of  the 
resistors  Ry  and  Rg,  at  the  last  reactor  node  point  inside  of  the  reactor, 
and  the  first  node  point  inside  of  the  reflector  need  not  be  adjusted.  How¬ 
ever,  if  the  interface  occurs  at  a  node  point  the  Ry  and  Rg  resistors  at 
this  point  must  be  adjusted.  To  do  this  the  respective  regions  are  considered 
to  overlap  at  that  node  point  and  the  following  value  is  used. 


Ry  (interface)  = 


2  Ry  (reactor)  x  R^,  (reflector) 
Ry  (reactor)  +  Ry  (reflector) 


(16) 


The  method  of  arriving  at  this  formula  is  shown  in  fig.  A-3.  If  the  inter¬ 
face  occurs  some  fraction  of  the  way  between  the  last  node  point  in  the 
reactor  and  the  first  node  point  in  the  reflector,  only  the  diffusion  resistor 
between  these  node  points  must  be  changed  if  R^  or  Rgfor  the  reactor  is 
not  equal  to  R^  or  R q  for  the  reflector.  In  the  case  where  the  interface  is 
at  one  half  of  the  distance  between  the  node  points  the  fast  network  diffusion 
resistor  are  equal  to 


R<£  (reactor)  plus  R$  (reflector) 
2  2 


and  the  thermal  flux  network  resistor  is  computed  in  the  same  way. 

Another  boundary  condition  which  must  be  simulated  is  the  insertion 
of  a  control  rod.  This  condition  may  require,  for  example,  that  the  thermal 
flux  be  zero  at  the  effective  radius  of  the  rod.  The  thermal  flux  is  made  to 
go  to  zero  by  grounding  the  0  network  at  points  in  that  network  which  corre¬ 
spond  to  the  dimensions  of  the  rod.  The  production  of  neutrons  is  zero  so 
again  no  current  feed  is  needed  in  that  section  of  the  0  network  which  was 
grounded.  In  this  way  any  number  of  rods  of  any  configuration  can  be  sim¬ 
ulated. 


G.  Conclusions 

* 

1.  A  network  can  be  constructed  so  as  to  simulate  the  reactor 
described  by  the  two-group  theory  equations. 


2.  The  network  is  one  of  lumped  constants  and  contains  a  finite 
distance  between  the  network  node  points.  These  finite  network  lengths 
have  a  definite  directional  relation  of  90  degrees  to  one  another. 

3.  The  accuracy  of  the  results  obtained  from  the  simulator  network 
depend  upon: 

a.  Selection  of  resistors  whose  values  closely  approximate  the 
calculated  values. 

b.  The  accuracy  with  which  the  current  feeding  equipment  fulfills 
the  theoretical  requirement:  I  =  AV g. 

c.  The  accuracy  of  instrument  reading. 

d.  The  error  involved  in  transforming  the  differential  equation  of 
the  reactor  into  finite  difference  form.  In  the  Taylor  expansion,  all  odd 
powered  terms  cancel  out  and  terms  of  the  fourth,  sixth,  etc.,  power  con 
stitute  the  error. 

The  total  error  resulting  from  all  of  these  effects  is  best  determined 
by  constructing  a  network  and  comparing  the  measured  results  with  those 
obtained  by  solving  the  reactor  equations  analytically. 


(a)  Square  Mesh  Network 
Fitted  To  Any  Boundary 


(b)  Details  of  One  Lattice  Unit  from 
Network  (a) 


PRO.  A't4oZ 


Generalized  Network 


Lattice  Unit  of  Reactor  Simulator 
Network 


P.P  P.  A-/4J3 


Fig. A- 


RRD.A'HO* 


28 


SECTION  B 


CONSTRUCTION  OF  EXPERIMENTAL  NETWORKS 


A.  A  one  dimensional  simulator  is  also  called  a  network  system. 

1.  General  Description:  The  first  test  network  constructed  by 

E.  J.  Wade  is  shown  on  fig.  B-5.  This  drawing  contains*the  complete  hook¬ 
up  including  the  network  and  the  current  feeding  amplifier  and  switch.  The 
network  is  the  equivalent  of  a  slab  reactor  in  two-group  theory  calculations. 
A  slab  reactor  is  one  in  which  the  neutron  flux  depends  on  only  one  coor¬ 
dinate.  An  8  /if  condenser  was  connected  between  each  point  of  the  0  net¬ 
work  and  ground.  The  condenser  together  with  the  rotating  switch  and  am¬ 
plifier  constitute  the  current  feeding  equipment  for  the  network  of  resistors. 
As  the  switch  revolves  it  connects  the  amplifier  input  and  output  succes¬ 
sively  to  corresponding  points  on  the  0  and  6  networks.  The  amplifier 
supplies  a  current  directly  proportional  to  the  input  voltage  and  independent 
of  the  load  resistance.  When  this  assembly  is  in  contact  with  any  pair  of 
network  points  the  current  charges  the  condenser.  The  condenser  then  con¬ 
tinues  to  feed  current  into  that  point  of  the  network  after  the  amplifier  and 
switch  assembly  moves  on  to  the  other  points.  This  process  is  referred  to 
as  “scanning.”  By  using  this  system  the  necessity  for  having  an  amplifier 
for  each  network  point  is  eliminated. 

The  main  purpose  of  this  network  was  to  study  the  feasibility  of  this 
scanning,  and  to  develop  a  suitable  amplifier. 

A  detailed  discussion  of  the  switch  and  the  amplifier  will  be  given 
in  Section  C. 

2.  Network  Resistance  Values:  (In  all  cases  the  available  standard 
resistor  size  closest  to  the  calculated  value  was  used.) 


reactor  length 
no.  of  lattice  units 


240  cm, 


=  10  cm 


A  standard  resistor  size  was  selected  for  R^  '=  2.2  megohms 


By  formula  10  Section  A,  with  T  -  250  sq  cm 

„  R^T  (2.2)  (250)  c  c  '  ,  , 

R  =  ,<£.  =  ^ - L  -  5.5  use  5.6  megohms. 

*  -> 


100 


A  standard  resistor  was  then  selected  for  Rq  =  0.047  megohms 
which  was  approximately  1  / 50  of  R^  then  by  formula  13a  of  Section  A 

-!*-  =  -%!=■=  46.8  =  k  Alth 
R#  .047  X 

Let  X^-£  =  Xtth*  Then  k  =  46.8  the  kappa  factor. 

With  L2  =  200,  by  formula  14a  of  Section  A 


R_  (1  -  L2  Xtth) 

Rg  _  L2  (  T  k  Xtf) 

0.047  _  100  (1  -  200  ) 

Rg  “  200  (  (250)  46.8) 

Rg  =  0.1  megohms 


For  the  reflector  let  T=  250  sq  cm,  the  same  as  in  the  reactor. 
Hence  the  Ky  in  the  reflector  =  the  R*y  in  the  reactor  when  the  lattice  unit 
length,  a,  and  R<£  are  not  changed. 

Let  L2  for  reflector  be  950  sq  cm 


0.047  100  (1  -  950  ) 

Rg  “  950  (  (250)  (46.8)  ) 

Rg  =  .47  megohms 

The  boundary  between  reactor  and  reflector  was  made  by  adjusting 
the  values  of  Rg  at  point  1  using  formula  (16).  See  fig.  C. 


R  (interface)  = 


2  (100,000)  (470,000) 
100,000  +  470,000 


=  164  000  ohms 


3.  Measurement  of  k;  The  X  network  system  is  operated  by  first 
starting  the  multicontact  switch.  The  output  current  of  the  amplifier  for 
zero  signal  is  not  exactly  zero  and  therefore,  a  small  initial  current  flows 
causing  voltages  to  be  established  at  the  various  points  of  the  0  network. 

The  small  initial  amplifier  current  corresponds  to  the  source  of 
neutrons  required  to  start  a  chain  reaction.  The  value  of  A  for  the  ampli¬ 
fier,  which  has  been  shown  to  represent  k  of  the  reactor  in  Section  A,  then 
feeds  a  current  representing  neutrons  into  the  network  as  dictated  by  this 
small  initial  source  producing  a  voltage  Vq  .  If  A  is  large  enough  to  rep¬ 
resent  super  critical  the  voltage  of  the  network  rises  just  as  the  flux  of  the 
reactor  increases  if  more  neutrons  areproduced  thanare  absorbed.  When 
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the  V$  approaches  a  convenient  value  such  as  100  volts  at  the  highest  point, 
A  is  reduced  until  ceases  to  rise.  This  indicates  that  just  enough  cur¬ 
rent  is  being  received  by  the  network  to  replace  the  current  lost  by  leakage 
to  ground.  This  represents  the  steady  state  condition  in  the  reactor  which 
calls  for  leakage  plus  absorption  of  neutrons  to  be  equal  to  neutron  produc¬ 
tion.  This  A  value  is  measured  in  amps  per  volt,  which  is  conductance,  and 


having  A, 


k  for  the  reactor  is  calculated  by  the  formula,  A  - — 

Kg 


+ 


k-l 

R'y 


for  the  close  approximation  due  to  “kappa  factor*5  since  the  large  value  of 
R  (5.6  megohms)  makes  the  last  term  small,  so  that  k  =  A.Rg- 


4.  Boundary  Conditions:  To  operate  the  network  as  a  core  without 
reflector  the  0  and  9  networks  were  grounded  at  the  network  points  corre¬ 
sponding  to  the  reactor  boundaries. 


To  simulate  the  effect  of  a  control  rod  the  Q  network  was  grounded 
at  the  point  in  that  network  at  which  the  flux  was  to  be  zero.  For  the  curve 
of  fig-  C-2  the  network  was  grounded  at  points  11L<,  12  and  11R.  Each  lat¬ 
tice  unit  represents  10  cm  of  core  material  hence  the  effect  of  a  1  x  1  x  20  cm 
rod  black  to  thermal  energy  neutrons  was  simulated.  In  fig®  B-l  and  B-2 
the  effective  rod  dimensions  are  1x1x0  cm  because  only  one  network 
point  was  grounded. 

5.  Experimental  Results:  In  the  test  runs  on  the  simulator  k  was 
measured  for  various  boundary  conditions.  The  k  was  then  calculated  by 
solution  of  the  reactor  equations  so  that  the  accuracy  of  the  measurement 
could  be  determined.  For  the  core  and  no  control  rods,  the  measured  k  was 
1.16;  the  calculated  k  was  1.078.  The  balance  of  the  results  were  taken 
with  the  insertion  of  control  rods  simulated.  Various  dimensions  of  rod 
were  simulated  and  in  each  case  the  rod  represents,  in  a  slab  reactor,  a 
slab  rod  of  1  cm  x  1  cm  x  some  desired  thickness.  The  thickness  was  var¬ 
ied  in  each  case. 


Figure  B-l  describes  a  core,  and  a  rod  of  zero  effective  thickness 
at  the  center.  The  measured  k  was  1.24;  the  calculated  k  was  1.270. 


Figure  B-2  describes  a  core  with  a  rod  of  20  cm  thickness  at  cen¬ 
ter.  The  measured  k  was  1.32;  the  calculated  k  was  1.360. 

Figure  B-3  describes  a  core  with  two,  zero  effective  thickness, 
rods  at  points  40  cm  out  from  center.  The  measured  k  was  1.695;  the 
calculated  k  was  1.825. 

Figure  B-4  describes  one  half  of  the  core  with  reflector  and  no  rods. 
The  measured  k  for  this  core  was  1.151,  the  calculated  k  was  1.06. 


Each  of  the  above  figures  shows  the  plot  of  the  flux  distribution,  as 
simulated  by  the  network  voltages,  as  a  solid  line.  For  comparison  with 
this  the  calculated  flux  distribution  from  an  analytical  solution  of  the  reac¬ 
tor  equations  was  plotted  as  a  dotted  line  for  B-l  and  B-2. 


B.  Two  Dimensional  Simulator 


1.  General  Description;  The  simulator  was  built  to  represent  a 
cylindrical  reactor  of  length  1  cm,  and  91.44  cm  over-all  radius,  having 
60.96  cm  radius  reactor.  For  a  two  dimensional  reactor  neutron  flow  is 
considered  to  be  in  two  directions  only,  that  is,  the  flow  under  considera¬ 
tion  is  confined  to  a  plane.  The  network  then  represents  a  sheet  of  reac¬ 
tor  material  1  cm  thick  having  the  above  dimensions.  The  flux  distribution 
in  such  a  reactor  is  symmetrical  around  the  center  and  therefore  the  network 
was  built  to  represent  only  a  half  circle  of  reactor  material.  The  plan  view 
of  this  network  is  shown  on  fig.  B-9.  To  construct  the  network  the  circu¬ 
lar  reactor  section  was  first  drawn  to  scale  and  then  divided  into  square 
sections  as  outlined  by  the  dotted  lines.  In  order  to  approximate  a  circu¬ 
lar  area  with  a  group  of  small  square  area  units,  the  size  of  the  units  must 
be  small.  In  this  case  two  sizes  of  square  area  units  were  used,  so  that 
the  effect  of  changing  network  lattice  size  could  be  studied.  Having  laid 
out  the  square  sections  a  network  node  point  was  placed  at  the  center  of 
each  area,  and  the  points  were  joined  by  resistors.  Each  large  unit  rep¬ 
resents  a  volume  of  reactor  material  15.24  sq  cm  and  1  cm  thick  and  each 
small  unit  represents  a  volume  7.62  sq.cm  and  1  cm  thick. 

The  drawing  shows  apian  view  of  the  network  representing  one  en¬ 
ergy  group  whereas  the  entire  network  consists  of  two  such  arrangements 
of  resistors  with  resistors  R-y  between  corresponding  node  points  of  the  0 
and  ©networks  and  resistors  Rg  to  ground  from  each  of  the  6  network  node 
points . 


2.  Current  Feeding  Equipment:  There  are  60  node  points  in  the 
section  of  network  which  represents  the  reactor.  Each  of  these  0  network 
node  points  must  receive  a  current  feed  proportional  to  the  0  network  volt¬ 
age  at  that  point  and  also  proportional  to  the  volume  of  reactor  material 
represented  by  the  lattice  unit  containing  that  point.  If  the  same  amplifier 
is  used  for  both  sized  volumes,  the  time  the  amplifier  feeds  each  point  must 
be  proportional  to  the  area  represented.  This  was  accomplished  by  feeding 
each  large  lattice  unit  four  times  as  often  as  the  small  units  were  fed.  The 
large  lattice  represents  four  times  the  volume  of  reactor  material  repre¬ 
sented  by  the  small  lattice. 

There  are  45  feed  points  for  small  volumes  and  15  feed  points  for 
large  volumes  so  a  switch,  having  (4  x  15)  plus  45  or  105  contact  pairs  was 
used.  It  was  decided  to  build  these  contacts  on  two  separate  disks,  each 


disk  having  70  pairs  of  contacts.  The  construction  of  this  switchis  covered 
in  Section  C.  It  is  immediately  clear  that  if  two  switches  are  used,  two 
amplifiers  will  be  necessary,  one  for  each  group  of  points  fed  by  a  switch. 
These  amplifiers  must  be  similar.  It  should  be  noted  that  even  though  the 
105  pairs  of  contacts  could  have  been  built  on  one  disk  in  this  particular 
case,  in  the  construction  of  three  dimensional  simulator  it  would  be  neces¬ 
sary  to  have  several  switches  and  amplifiers.  Two-disk  construction  was 
used  so  that  the  problem  of  matching  the  characteristic  curves  of  more  than 
one  amplifier  could  be  studied. 

3.  Calculation  of  Resistor  Values:  It  was  pointed  out  that  the  lat¬ 
tice  size  is  directly  related  to  the  volume  of  reactor  represented  by  that 
network  section.  It  is  clear,  then,  that  if  a  value  of  R0  is  selected  and 
used  throughout  the  network  the  core  volume  represented  by  a  network  lat¬ 
tice  unit  is  determined  by  the  value  of  R-y  at  the  node  point  of  that  section 
of  lattice.  Let  the  pile  constants  for  the  reactor  be  T  =  258,  L2  =  121  and 
Xtf  =  Xtth.  The  value  of  R0  was  first  selected  to  be  2.25  megohms  and 
also  the  value  of  R@  ,  0.047  megohms. 

Selecting  the  values  of  R0  and  R@  for  the  small  lattice  unit  fixed 
the  value  of  R0  and  RQ  for  the  larger  lattice  units.  When  changing  the  lat¬ 
tice  unit  size  the  diffusion  resistance  per  unit  area  must  be  kept  the  same. 
Referring  to  fig.  B-9,  for  example  at  point  11,  the  resistance  in  either  the 
x  or  y  directions  of  a  small  unit  is  R  for  the  area  h  (7 .62  x  7.62  sq  cm) 
enclosed  by  the  dotted  line  around  point  11.  Compare  this  area  with  the 
four  points  11,  21,  22,  12  which  form  a  square  and  are  together  enclosed 
by  dotted  lines  of  2h  cm  on  a  side.  The  area  contained  by  these  lines  is 
therefore  4h2.  These  dotted  lines  always  bisect  a  resistor,  so  in  the  x 
direction  starting  at  a  left  hand  dotted  area  it  is  found  that  there  is 
R/2  +  R  +  r/z  or  2R  ohms  along  the  line  of  21  and  11  and  the  same  2R  ohms 
along  the  line  of  22  and  12.  The  parallel  combination  is  R,  the  resistance 
through  which  current  must  pass  in  going  from  the  left  hand  boundary  of 
the  four  units  to  the  right  hand  boundary  in  the  direction  x.  The  same  can 
be  said  of  the  direction  y  at  90  degrees  to  x.  It  follows  logically  then,  that 
for  the  area  of  4h2  a  single  resistance  R  could  be  installed  between  the  two 
boundaries  in  either  the  x  or  y  direction  and  the  current  passing  across 
the  section  would  not  be  changed.  This  is  what  was  done  in  the  large  lat¬ 
tice  units.  For  example,  at  point  022  r/2  +  R/2  =  R  for  the  resistance  R 
in  either  the  x  or  y  directions.  Only  the  x  and  y  components  of  current 
diffusion  which  has  been  demonstrated  in  Section  A  in  the  development  of 
the  network  formulas  are  considered.  Current  which  diffuses  in  any  inter¬ 
mediate  direction  between  x  and  y  is  broken  up  into  its  x  and  y  components. 
Consider  that  there  exists  a  voltage  at  point  022  greater  than  that  at  044. 

A  current  cannot  flow  directly  from  point  022  to  point  044  if  a  direct  con¬ 
ductor  connection  does  not  exist  between  them.  However,  a  current  does 
flow  between  these  points  through  the  balance  of  the  network,  because,  by 
the  superposition  network  the'orem,  the  flow  of  current  in  any  network  due 


to  any  one  voltage  is  independent  of  other  voltages  existing  on  the  network. 
It  is  important  therefore  that  in  designing  a  network  the  X  and  Y  direction 
resistances  per  unit  square  area  be  maintained  the  same,  no  matter  what 
size  area  is  considered.  Since  the  distances  in  the  X  and  Y  directions  are 
equal  for  square  mesh,  a  =  b  =  c=  h  =  7.62  cm  for  the  small  lattice  units, 
from  10 

R<2>  h2 
Ry  T 


R 


7  " 


(2.25)  .  (258) 
(7.62)2 


=  10  megohms 


Here  again  as  in  the  slab  reactor  simulator,  a  M  kappa  factor”  is 
introduced  so  that  the  0  network  voltage  will  be  small. 

By  formula  13a 


R^  _  2.25  _  _  K  ^  th 

R0  .047  Xtf 

Then  by  formula  14a 

Re  _  p,  Rz  x  tf  1 

Rg  _  L2  L  “  Tk  X  tthj 


121  x  .047 _ 

Kg  ~  (7.62)2  x  (1  -  121  i _ ) 

258  '  47.87) 


=  0.09795  megohms 
use  0.1  megohms 


For  the  large  lattice  units  h  =  15.24  cm 


F  rom  10 

R<£  _  h* 

Ry  T 

Since  T  does  not  change,  and  R  is  the  same  for  both  sized  lattice 

units, 


(large  mesh)  _  h2  (small  mesh)  (7.62)2  _  1 

Ry  (small  mesh)  h2  (large  mesh)  (15. 24)2  4 


R,y  (large  mesh) 


R-v  (small  mesh) 

— I —  —  *  - .  ■  -  6,0  megohms 
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Similarly 

Rg  (large  mesh)  _  Rg;( small  mesh)  =  25>000  ohms 

4 

Here  the  nearest  available  size  was  24,000  ohms  which  is  the  value 
used  in  the  network. 

For  the  reflector,  letting  T  =  150  sq  cm,  L.2  =  100  sq  cm  and  using 
the  diffusion  resistor  values  of  R0  and  R Q  which  are  now  fixed  values  the 
following  is  obtained. 


R_.  =  (2-25)  (15°).  _  5.813  megohms  for  the  small  mesh  for  which  the 
1  (7.62)2 

value  5.6  megohms  was  used. 


5.6 

R  (large  mesh)  =  — 
1.5  megohms  was  used. 


=  1.4  megohms  for  which  the  standard  value 


„  ,  (1000)  (0.047) 

Rg  (small  mesh)  =  j~z  _ 


(7.62)2  (1  -  1000  .  _i 


=  1  megohm 


150  47.87 


Rg  (large  mesh)  =  0.24  megohms 

4.  Boundary  Conditions:  On  fig.  B-9  it  can  be  seen  that  the  bound¬ 
ary  between  the  reflector  and  reactor,  represented  by  the  heavy  dotted  line, 
passes  between  the  node  points  of  the  network.  As  a  result  of  this  choice 
it  is  not  necessary  to  adjust  the  value  of  R-y  or  Rg  at  this  boundary  as  was 
necessary  in  the  slab  reactor  simulator  when  the  boundary  passed  through 
a  node  point.  The  boundary  of  the  reflector  is  grounded  for  both  the  0  and 
0  networks. 

The  method  of  making  the  connections  between  the  large  and  small 
mesh  sections  must  still  be  justified,  (refer  to  fig.  B-9).  In  the  x  or  y 
direction  the  value  of  R  for  the  small  lattice  corresponds  to  the  distance  h 
whereas  in  the  large  lattice  it  corresponds  to  2h.  Hence,  the  distance  h 
in  the  large  lattice  requires  the  value -^v  In  the  small  lattice  the  distance 
^  corresponds  to  Examining  the  adjacent  points  0,  1  and  01  on  the  ^ 
fig.  B-9  it  can  be  seen  that  the  distance  from  0  or  1  to  01  is  the  sum  of  y 
for  the  large  mesh  plus  ^  for  the  small  lattice.  If  the  resistor  j-  for  the 
large  lattice  is  split  into  two  resistors  in  parallel,  each  of  value  R,  and 
these  resistors  each  placed  in  series  with  the  -jr  required  for  each  of  the 
small  lattice  points,  the  total  resistance  between  either  0  or  1  and  01  is 
1.5  R.  This  is  the  value  used.  This  connection  corresponds  to  the  physical 
*  concept  of  neutron  diffusion  inasmuch  as  the  neutrons  at  01  which  diffuse 


* 


in  the  X  direction  must  reach  0  and  1.  The  reverse  is  also  true  because 
neutrons  from  points  0  and  1  must  flow  to  01. 

5.  Control  Rods  in  Reactor:  The  network  as  designed  will  allow 
simulating  the  insertion,  into  the  core,  of  two  sizes  of  rods,  black  to  ther¬ 
mal  neutrons,  one  7.62  sq  cm  and  the  other  15.24  sq  cm  by  grounding  the  0 
network  at  the  boundaries  of  a  0  network  lattice  unit. 

If,  for  example,  a  rod  7.62  sq  cm  is  desired  at  point  10,  the  resis¬ 
tors  R  are  disconnected  from  between  points  10  and  20,  10  and  11,  10  and  0. 

■p 

Resistors  of  value  ^  ohms  are  then  connected  to  ground  from  each  of  points 
20,  11  and  0.  The  current  feed  to  point  1  0  in  the  0  network  must  then  be 
zero  to  simulate  no  production  of  neutrons  in  the  reactor  volume  represented 
by  the  dotted  square  around  10.  Likewise  point  022  can  be  treated  in  the 
same  manner.  By  placing  resistors  of  ^  ohms  between  02,  024,  042,  020 
and  ground,  a  rod  15.42  sq  cm  is  simulated. 

The  location  of  these  rods  is  confined  to  areas  around  existing  node 
points  of  the  network.  Generally,  however,  a  method  is  desired  of  simu¬ 
lating  a  rod  of  any  given  size  at  any  point  in  the  reactor.  To  do  this  there 
must  exist  at  that  point  a  network  lattice  unit  which  represents  the  cross 
sectional  area  of  rod  to  be  simulated.  The  method  of  changing  network  lat¬ 
tice  size  described  above  can  be  extended  through  subsequent  reductions 
until  the  size  of  area  simulated  is  very  nearly  the  area  of  control  rod  de¬ 
sired.  The  subsequent  reductions  in  area  size  are,  starting  with  an  area 
h2 ;  h^  h2  h2  etc.  The  value  of  R  and  Rg  for  each  of  these  points  must.be 
4  3  16’  64 

adjusted  in  accordance  with  the  formulas  in  Section  A  for  each  new  section 
length,  equal  to  etc.  At  the  same  time  the  current  feeds  to  the  0 

network  must  be  changed  so  that  the  current  fed  to  each  point  is  propor¬ 
tional  to  the  area  represented.  This  current  feed  provision  is  the  most 
difficult  to  satisfy.  There  are  various  methods  by  which  this  can  be  done 
and  they  are  covered  in  Section  D. 

6.  Experimental  Results:  All  experimental  work  to  date  has  been 
done  with  the  network  simulating  a  reactor  without  control  rods.  For  this 
type  of  reactor  the  calculation  of  the  flux  distribution  and  k  by  analytical 
solution  of  the  reactor  equation,  is  not  a  very  difficult  problem. 

The  process  of  operating  the  two  dimensional  simulator  involves 
inserting  a  voltage  source  at  some  point  on  the  0  network  and  thereby  set¬ 
ting  up  some  distribution  of  voltage  in  the  network.  This  simulates  the 
presence  of  a  neutron  source.  The  switch  and  amplifier  are  then  started 
and  because  a  current  I  =  A Vq  is  then  fed  into  the  0  network  points,  the 
voltage  distribution  becomes  a  function  of  the  network  resistor  values  and 
the  value  of  amplifier  transconductance  A.  In  the  simulator  used  for  the 
experiments  two  similar  amplifiers  were  used  to  feed  the  network  through 


the  two  separate  switches.  The  voltage  distribution  is  seriously  affected 
by  any  dissimilarity  in  the  value  of  A  for  each  of  these  amplifiers.  The 
characteristic  input  voltage -output  current  curves  for  these  amplifiers 
must  coincide  throughout  the  range  of  current  feed.  Similarly,  the  con¬ 
tact  time  for  each  switch  must  be  the  same  because  if,  for  some  reason, 
one  switch  fed  current  for  somewhat  longer  intervals  than  the  other,  that 
network  section  would  receive  more  current.  This  would  result  in  an 
improper  voltage  distribution  on  the  network,  because  the  current  feed  (i) 
would  not  even  approximate  the  desired  or  design  value. 

The  two  amplifiers  displayed  a  tendency  to  have  a  shift  in  their 
zero  signal-current  output  and  the  result  was  a  change  in  the  slope  A  of 
their  characteristic  curves.  For  this  reason  it  was  found  that  adjusting  the 
amplifiers  to  coincide  before  a  run  was  not  sufficient  and  that  a  continual 
check  of  this  transconductance  was  needed.  In  practice  this  check  consisted 
of  installing  a  milliammeter  and  a  voltmeter  at  representative  network  feed 
points  in  each  section,  fed  by  an  amplifier.  The  current  per  volt  at  these 
points  was  measured  and  compared  by  the  operator  so  that  he  could  adjust 
the  A  of  the  amplifiers  to  coincide  at  one  point  on  their  characteristic 
curves.  This  operation  is  covered  more  fully  in  Section  C-5.  Inaccuracies 
due  to  the  switch  contact  time  could  not  be  adjusted  in  all  cases. 

The  experimental  results  are  therefore  somewhat  dependent  on  the 
operator's  observations  during  a  run.  It  will  be  necessary  in  future  work 
to  eliminate  this  accuracy  factor  by  developing  a  more  stable  amplifier  and 
a  more  accurate  switch. 

The  results  obtained  from  the  simulator  are  plotted  on  fig.  B-6  and 
B-7.  The  flux  distribution  obtained  by  an  analytical  solution  of  the  reactor 
equation  is  plotted  on  the  same  graph  for  comparison.  The  k  measured 
using  the  simulator  was  1.263  and  the  analytical  value  was  1.266. 
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PLOT  OF  4>f  FOR  SIMULATOR  NETWORK  FED  BY  AMPLIFIER  AND 
SCANNING  SWITCH.  LARGE  LATTICE  SECTION. 
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FIG.  B-6  PPD-A-1397 


PLOT  OF  £ f  FOR  SIMULATOR  NETWORK  FED  BY  AMPLIFIER  AND 
SCANNING  SWITCH.  SMALL  LATTICE  SECTION. 
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FIG.  B-7  PPD-A-1398 


SECTION  C 


CURRENT  FEEDING  EQUIPMENT 


A.  General 

The  condensers,  scanning  switch  and  amplifier  constitute  the  current 
feeding  equipment  for  the  resistor  network.  The  operation  of  this  assembly 
can  be  studied  by  resolving  the  system  into  the  simple  circuit  of  fig.  C-l 
(a). 

G  is  a  constant  current  generator  which  represents  the  amplifier. 

The  amplifier  has  virtually  a  constant  current  output  for  any  particular  net¬ 
work  point  when  the  network  voltages  and  Vg  are  nearly  constant. 

C  is  the  10  /if  condenser  between  the  0  network  node  points  and 

ground. 

R  is  the  resultant  network  resistance  from  any  node  point  to  ground 
for  a  steady  voltage  distribution.  This  resistance  changes  as  the  network 
voltages  change  since  the  effective  resistance  to  current  flow  at  any  par¬ 
ticular  point  is  a  function  of  not  only  R-y  but  also  of  the  voltage  difference 
between  the  node  point  being  fed  and  the  neighboring  node  points  of  the  net¬ 
work. 

S  is  one  pole  of  the  scanning  switch  which  has  an  open  time  of  0.257 
sec  (to)  and  a  closed  time  of  0.004  sec  (tc)  for  the  small  lattice  units  of  the 
two  dimensional  network. 

I  is  the  current  which  flows  when  the  switch  is  closed  so  that 
I  =  ix  +  i2  with  R/v'7  megohms  ii  »*2 

i  is  the  loop  currents  which  flows  when  the  switch  is  open.  The 
current  flow  at  point  A  is  plotted  in  fig.  C-l  (b).  If  I  is  constant  the  average 
value  of  V 


Vare  = 


R  I  tc 
tc  +  to 


The  voltage  rises  to  a  maximum  at  tl  then  falls  from  tx  to  t2  and  rises  again 
from  t2  to  t3  etc. 


With  RC  large  the  voltage  is  nearly  constant  and  in  practice  it 
was  found  that  the  variation  in  voltage  due  to  scanning  was  0.25  volts,  at 
100  volts  on  the  0  network,  or  a  variation  of  .25  per  cent.  This  is  well 
within  the  limits  of  experimental  accuracy. 
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At  steady  state  i  is  nearly  equal  to  i2  and 

.  _  Xlk-  —  1  d  1 

1  "  R  C  C 


With  I  constant  over  the  period  o  to  tx 


i  =  -  fl 

c  Jo 


dt  = 


I  t  c 
c 


Thus  the  current  i  is  seen  to  be  a  direct  function  of  the  time  through¬ 
out  which  the  switch  is  closed.  If  this  time  is  different  on  subsequent  switch 
closings  i  and  (as  a  result)  V0  will  change.  Furthermore,  if  the  time  tc 
were  different  for  the  different  contacts  of  a  scanning  switch  the  current 
fed  to  the  respective  network  feed  points  would  not  be  a  function  of  the 
amplifier  current  alone  but  would  vary  with  the  contact  time.  It  is  therefore 
important  that  the  scanning  switch  contact  time  be  constant  so  that  the 
amplifier  can  perform  its  function  properly. 

For  proper  scanning  of  the  large  lattice  unit  section  of  the  network, 
the  intervals  between  feeds  must  be  as  uniform  as  possible  or  an  additional 
voltage  variation  will  result.  To  understand  this,  consider  the  action  of 
scanning  a  network  point  four  times  in  one  revolution  of  the  switch.  In  fig.  C  —  l 
(c)  the  time  base  is  one  revolution,  tc  and  to  are  the  closed  and  open  times 
respectively  for  a  contact,  and  1,  Z,  3,  4  are  subscripts  denoting  the  four 
separate  contacts  which  feed  current  to  the  network  point. 

The  voltage  rise  due  to  the  current  feed  of  contact  number  1  is  0.Z5 
volt.  For  the  small  lattice  units  the  voltage  decay  is  very  slow  while  for  the 
large  lattice  units  the  decay  is  four  times  as  fast.  In  fig.  C-l  (c)  the  four 
feed  contacts  have  been  assumed  to  be  spaced  evenly,  and  therefore  the 
peak  of  the  saw  tooth  wave  is  always  at  the  same  voltage.  In  the  two  di¬ 
mensional  simulator  the  contacts  were  actually  spaced  as  shown  on  fig.  C-Z. 
Contacts  numbers  1  to  Z,  Z  to  3,  3  to  4  have  13  segments  between  each,  and 
between  numbers  4  and  1  there  were  31  segments.  The  result  of  this  unequal 
spacing  was  a  variation  from  maximum  to  minimum  voltage  of  about  0.35 
volts.  In  future  construction  more  evenly  spaced  contacts  would  be  desirable. 


B.  Scanning  Sequence 

The  sequence  in  which  the  feed  points  of  the  networks  are  fed  by 
the  amplifier  has  a  significant  effect  on  the  voltage  distribution.  It  was 
found  in  the  work  with  the  slab  reactor  simulator  that  scanning  the  points 
in  rotation  from  one  end  of  the  network  to  the  other  resulted  in  an  improper 
voltage  distribution  whereas,  if  the  network  was  fed  as  shown  on  fig.  B-5, 
the  voltage  distribution  more  nearly  simulated  the  calculated  flux  distribution. 


Feeding  adjacent  network  node  points  consecutively  is  not  desirable 
because  the  network  voltage  of  points  abound  the  point  which  has  just  been 
fed  are  at  the  peak  of  the  saw  tooth  of  fig.  C-l.  Since  the  current  fed  to  a 
point  is  proportional  to  the  voltage  at  that  point,  the  current  fed  will  be 
higher  than  it  would  have  been  had  some  time  elapsed  for  the  voltage  to 
decay  before  the  current  feed  was  applied. 

In  the  two  dimensional  simulator  the  scanning  sequence  was  selected 
so  that  the  points  were  fed  almost  at  random.  The  sequence  is  shown  in 
fig.  C-3  for  the  large  lattice  section.  The  sequence  for  the  section  con¬ 
taining  small  mesh  is  as  follows: 

Switch  segment  number  1  was  connected  to  network  point  number  73: 
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C.  Time  Factors 

The  time  required  to  charge  the  condensers  connected  to  the  network 
determines  the  time  required  to  stabilize  the  network  at  any  desired  voltage. 
This  time  factor  has  no  connection  with  the  reactor  transients. 
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When  the  network  was  completely  de-energized  and  all  the  capacitor  s 
were  discharged  it  was  found  that  the  time  required  to  establish  steady  state 
network  voltage  of  100  volts  peak  value  was  about  15  minutes.  This  15 
minutes  was  not  a  fixed  figure  but  the  average  time  required  by  the  operators 
who  worked  with  the  simulator.  If  the  amplifier  were  applied  to  a  completely 
de-energized  network,  where  there  was  no  Vq  voltage,  the  theoretical  cur¬ 
rent  fed  is  zero,  just  as  in  a  reactor,  no  fissions  take  place  unless  a  neutron 
source  is  present.  Actually,  the  current  output  of  the  amplifier  was  never 
exactly  zero  so  some  current  would  flow,  in  the  order  of  5x10  amps, 
and  the  cycle  of  production  could  start.  The  voltages  increase  as  long  as 
A  is  above  the  critical  value,  in  which  case  a  super  critical  k  exists  for  the 
reactor  simulated.  A  must,  of  course,  be  reduced  to  the  critical  value  when 
the  desired  peak  voltage  has  been  established.  If  this  small  initial  current 
were  depended  upon  in  establishing  the  network  voltages  for  a  run  the  time 
consumed  in  reaching  a  peak  voltage  of  100  volts  would  be  very  long  unless 
A  was  set  far  above  the  critical  value.  With  A  set  only  slightly  super  crit¬ 
ical  a  long  time  is  required  and  the  tendency  is  for  the  operator  to  set  A 
very  high  to  speed  up  the  charging  process.  If  this  is  done  the  voltages 
rise  rapidly  as  the  Vq  voltages  become  larger.  To  prevent  the  Vq  voltages 
reaching  the  maximum  voltage  output  of  the  amplifier  A  must  be  reduced 
to  slightly  above  critical  as  soon  as  the  voltages  begin  to  rise  rapidly. 

The  amplifier  characteristic  is  not  I  =  A  Vq  at  load  voltages  above 
about  180  volts  and  the  current  feed  requirement  is  not  being  met  by  the 
amplifier  in  this  region  of  saturation  so  the  system  does  not  simulate  the 
reactor.  It  was  found  that  the  best  way  to  start  up  the  simulator  was  to 
apply  a  Vq  voltage  at  the  network  node  point  where  the  peak  voltage  would 
logically  occur  and  hold  the  value  of’  A  just  slightly  above  critical  during 
the  buildup  time. 

The  possibility  exists  that  the  Y 0  voltages  will  rise  to  the  voltage 
saturation  point  of  the  amplifier  during  a  run  as  well  as  when  starting  up. 
With  the  network  energized  and  A  set  at  the  critical  value  the  system  is 
theoretically  stable.  Actually  this  is  not  the  case.  A  slight  change  in  A 
due  to  a  shift  in  amplifier  characteristics  caused  the  network  voltage  to 
fall  or  rise.  If  the  value  of  A  was  returned  to  the  critical  value  immediately 
the  voltages  fell  to  zero  or  rose  to  the  maximum  voltage  dictated  by  the 
amplifier.  This  rise  or  fall  was  slow  if  only  a  small  shift  from  critical 
took  place  but  none  the  less  the  voltages  changed  and  the  voltage  distribution 
could  not  be  read  accurately.  To  prevent  this  voltage  drift  it  was  imperative 
therefore  to  include  in  the  amplifier  a  means  of  providing  a  small  step 
change  in  A  which  could  be  controlled  by  the  network  voltage.  If  the  Vq 
voltage  rose  above  the  desired  peak  value  a  set  of  voltage  sensitive  relay 
contacts  were  opened  and  a  resistance  was  inserted  in  the  A  control  circuit. 
With  a  fall  in  voltage  the  reverse  took  place.  The  critical  A  value  then  was 
between  the  high  and  low  value  applied  through  this  control.  The  value  of 
A  must  be  determined  with  instrument  readings  taken  at  the  network, 


namely,  average  I  fed  to  a  point  divided  by  the  average  Vq  for  that  same 
point.  The  total  change  in  I  and  Vq  for  the  network  used  in  the  tests  was 
so  small  that  hardly  any  change  was  noticed  on  the  instruments.  It  was 
therefore  not  difficult  to  take  readings  which  represent  the  actual  current 
per  volt  fed  to  the  network. 

The  operator  must  have  the  ability  to  adjust  the  network  to  operate 
at  a  desired  peak  voltage,  hold  it  there  by  adjusting  the  voltage  sensitive 
control  relay,  and  then  read  the  instruments  properly.  This  ability  is  not 
difficult  to  acquire  through  a  reasonable  amount  of  experience  with  the 
simulator  system. 


D.  Scanning  Switch 

The  original  scanning  switch  for  the  slab  reactor  simulator  was 
made  by  inserting  pins  in  a  lucite  sheet  about  l/z  in.  thick.  Twenty -four 
pins  were  spaced  evenly  around  a  circle  and  a  contact,  carried  on  a  motor- 
driven  arm,  was  pivoted  at  the  center  of  the  circle.  See  fig.  C-4. 

The  contact  spring  was  set  so  that  it  did  not  touch  two  adjacent  pins 
at  the  same  time. 

It  has  been  pointed  out  that  uniformity  of  contact  time  is  very  im¬ 
portant.  Several  factors  contributed  to  the  poor  uniformity  of  this  switch: 
the  contact  bounced  when  the  spring  hit  a  pin;  the  spring  tension  changed 
during  operation;  all  the  pins  were  not  set  absolutely  vertical  to  the  mount¬ 
ing  sheet,  and  were  not  evenly  spaced.  The  0  network  connected  to  one 
ring  of  pins  and  the  Q  network  to  the  other  as  shown  on  fig.  B-5, 

An  improved  design  was  developed  and,  since  the  two  dimensional 
simulator  was  in  view  at  the  time,  it  was  made  with  70  pairs  of  contacts 
on  each  of  two  contact  discs.  The  complete  assembly  is  shown  on  fig.  C-5, 
which  is  a  photograph  of  the  switch  and  network  assembly.  Figure  C-6 
shows  the  detail  of  the  contact  assembly.  A  machine  screw  held  each 
contact  in  place  on  the  lucite  disc.  A  wire  was  connected  from  each  contact 
to  the  network.  The  outer  row  of  contact  was  connected  to  the  0  network 
points  and  the  inner  row  was  connected  to  corresponding  0  network  points. 
The  lucite  between  each  contact  was  machined  flush  with  the  metallic 
contact  surface.  A  revolving  contact  arm  carried  two  spring  steel  contacts 
which  slid  over  the  rings,  contacting  a  corresponding  pair  simultaneously. 
Electrical  connection  was  made  to  these  moving  contacts  through  two  slip 
rings  which  were  in  turn  connected  to  the  amplifier  input  and  output.  One 
side  of  the  network  and  amplifier  was  connected  to  ground  to  complete  the 
circuit. 
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The  contact  time  uniformity  for  this  switch  was  better  than  for  the 
original  switch  but  was  still  far  from  the  uniformity  required.  Figures  C-7 
and  C-8  are  photographs  of  the  switch  contact  operation  as  shown  on  a 
cathode  ray  oscillograph.  In  these  pictures  a  slow  vertical  beam  sweep 
was  coupled  with  a  fast  horizontal  sweep  to  make  successive  horizontal 
lines,  which  start  at  the  bottom  of  the  scope  screen  and  move  upward. 

The  switch  contacts  when  closed  caused  the  blank  spaces  in  these  lines. 

The  break  in  the  bottom  line  represents  one  contact  and  the  next  break 
above  it  is  the  next  contact  on  the  disc.  The  switch  was  turning  at  about 
230  rpm  for  fig.  C-7  and  about  115  rpm  for  fig.  C-8. 

The  large  variation  in  contact  time  is  easily  seen  from  these  two 
pictures.  Although  they  do  not  contain  all  the  contacts  they  show  a  repre¬ 
sentative  sample  of  the  contact  action.  The  contact  time  variation  was 
over  100  per  cent  from  the  shortest  to  the  longest  contact  time  shown  in 
the  photograph.  Poor  machine  shop  accuracy  caused  this  variation.  It  is 
not  likely  that  this  type  of  switch  can  be  constructed  accurately  enough  to 
meet  the  needs  of  the  simulator. 

Several  ideas  were  advanced  for  further  improvements  and  some 
test  work  has  been  done  with  a  view  to  developing  a  better  switch.  This 
work  is  described  in  Section  D  of  this  report. 


E.  Amplifier  ' 

It  has  been  pointed  out  that  a  current  must  be  introduced  at  each 
network  feed  point  in  the  0  network  which  is  proportional  to  the  voltage  of 
the  corresponding  Q  network  point.  The  same  amplifier  was  used  for  both 
the  slab  reactor  and  the  two  dimensional  simulator. 

The  ideal  input-output  relation  for  this  amplifier  is  a  straight  line 
I  (out)  =  A  Vq  (signal). 

*  * 

The  circuit  developed  by  E.  J.  Wade  is  shown  on  fig.  C-9-  This 
drawing  covers  the  amplifier  circuit  alone.  The  equipment  included  a 
standard  d-c  power  supply  which  is  not  shown.  An  additional  circuit  called 
a  compensator  was  built  onto  the  same  chassis  but  never  used  in  operation 
on  the  simulator.  This  circuit  is  not  shown  in  this  report.  The  purpose  of 
this  circuit  was  to  compensate  for  the  zero  drift  of  the  amplifier  by  chang¬ 
ing  the  voltage  across  R5  of  the  amplifier  input  circuit.  Compensation  for 
zero  drift  in  the  amplifier  was  needed  but  this  particular  method  did  not 
prove  satisfactory. 

The  power  supply  was  not  as  ripple  free  as  was  desirable  and  the 
circuit  of  fig.  C-9  was  changed  to  that  shown  on  fig.  C-10,  in  which  the 
three  OB-2  voltage  regulator  tubes  were  installed  in  place  of  the  voltage 
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divider.  In  addition  the  500,  0.5  megohm  potentiometer  was  inserted  in 
place  of  the  fixed  resistor  in  fig.  C-9.  This  gave  more  range  to  zero 
signal  adjustment.  For  operation  on  the  network  the  amplifier  was  so 
adjusted  that,  at  zero  voltage  input,  the  output  current  was  approximately 
5  x  10“6  amp  for  all  settings  of  the  A  control.  The  desired  zero  signal 
current  was  zero,  but,  since  this  point  of  the  characteristic  curve  was 
not  used  in  operation,  the  small  zero  current  was  tolerated.  In  setting 
the  zero  current  it  was  important  not  to  suppress  the  zero  point.  If  this 
was  done  a  considerable  voltage  signal  was  required  to  produce  any  cur¬ 
rent  flow.  Zero  must  be  set  where  the  output  is  just  5  x  10~6  amp  at  zero 
voltage  signal.  The  test  curves  of  this  amplifier  have  been  plotted  on 
fig.  C  —  1 1 .  From  these  curves  it  can  be  seen  that  the  effect  of  changing 
load  is  small  but  that  the  current  flow  for  a  given  signal  voltage  is  not 
entirely  independent  of  the  load.  For  signal  voltages  above  0.1  ^olts  the 
curve  was  almost  a  straight  line.  The  maximum  voltage  output  for  the 
amplifier  was  197  volts  and  it  was  necessary,  therefore,  to  operate  the 
network  at  voltages  well  below  this  value.  The  amplifiers  were  usually 
adjusted  with  the  output  feeding  current  to  ground.  This  is  nearly  the  cor¬ 
rect  loading  because  the  10  juf  capacitors  connected  to  each  node  point 
make  the  effective  load  on  the  amplifier  almost  purely  capacitive.  The 
main  difficulty  with  this  amplifier  was  the  fact  that  the  zero  setting  shifted 
and  would  not  remain  constant  long  enough  to  complete  a  run.  With  the 
zero  set  at  5  x  10~6  amp  at  the  start  of  a  run  it  was  found  to  be  50  x  10“6 
after  a  run.  Shifting  the  zero  point  caused  the  characteristic  curve  to 
shift.  This  affected  the  stability  of  the  network  voltages.  The  net  effect 
of  this  zero  current  change  on  the  accuracy  of  the  simulator  could  not  be 
evaluated  during  the  experimental  runs  because  of  the  difficulties  which 
existed  in  the  scanning  switch.  An  amplifier  is  required  which  minimizes 
the  effect  of  changes  in  tube  characteristics.  Such  an  amplifier  would  give 
more  stable  simulator  operation.  From  the  foregoing  discussion  of  the 
characteristic  of  one  amplifier  it  can  be  easily  understood  that  some  dif¬ 
ficulty  would  be  encountered  in  trying  to  use  two  amplifiers  to  feed  different 
sections  of  the  same  network.  For  such  an  operation  the  characteristic 
curves  for  the  two  amplifiers  should  coincide  for  all  values  of  the  A  control 
setting.  Here  again  any  change  in  zero  current  caused  the  curves  not  to 
coincide  and  the  network  section  fed  by  each  amplifier  did  not  receive  the 
same  current  as  was  specified  in  the  simulator  theory,  (i  =  A  Vq  ). 

Suppose  the  characteristic  curves  coincide  for  both  amplifiers  at 
start  of  a  run  and  are  shown  by  line  1  of  fig.  C-12.  Line  2  shows  the  effect 
of  a  change  in  zero  current  for  one  amplifier.  The  line  is  almost  parallel 
to  the  original  line  but  this  amplifier  then  supplies  a  current  per  volt  which 
is  higher  at  all  points  throughout  the  operating  range  than  its  mate.  Too  high 
a  voltage  is  then  developed  on  the  network  section  fed  by  this  amplifier  and 
the  voltage  distribution  on  the  whole  network  becomes  distorted,  the  peak 
voltage  being  shifted  toward  the  section  fed  by  the  amplifier  having  the  higher 
current  output  per  volt.  If  the  A  control  is  then  changed  for  this  amplifier, 
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the  slope  of  its  curve  will  change  as  shown  by  line  3,  causing  curves  for 
the  two  amplifiers  to  cross  or  coincide  at  one  point.  If  this  point  (X  on  the 
figure)  is  near  the  center  of  the  operating  range  the  characteristic  curves 
will  almost  coincide  over  the  operating  range  for  the  two  amplifiers.  In 
order  to  make  this  adjustment  during  a  run  an  ammeter  was  installed  in 
the  line  feeding  current  to  the  network  node  point  022  in  the  large  lattice 
section  and  22  in  the  small  lattice  section  (see  fig,  B-9).  The  current  per 
volt  for  these  two  network  node  points  was  then  adjusted  by  the  operator 
so  that  they  were  the  same. 


F.  Conclusions 


In  this  section  the  inherent  inaccuracies  of  the  current  feeding 
equipment  have  been  pointed  out.  No  evaluation  of  the  effect  of  each  was 
made  since  the  net  effect  of  all  errors  introduced  by  the  various  elements 
was  the  only  real  concern,  A  complete  investigation  of  all  possible  current 
feeding  methods  would  include: 

1.  Scanning  with  an  electronic  multivibrator  circuit  in  place  of  the 
mechanical  switch  with  an  amplifier  suitable  for  that  type  of  scanning. 

2.  Feeding  current  with  manually  adjusted  variable  resistors  con¬ 
nected  to  each  feed  point  from  a  source  of  voltage.  This  is  the  Iteration 
Method  described  in  Section  D, 

3.  Development  of  an  automatic  system  of  adjusting  the  variable 
resistors  (similar  to  2)  using  a  servo  system. 

4.  A  vacuum  tube,  for  each  feed  point,  whose  conductance  is  changed 
by  scanning  the  grid  circuit  of  the  tube  with  a  voltage  signal  amplifier . 

5.  Combinations  of  the  above  methods.  For  example  the  amplifier 
and  scanning  switch  method  can  be  used  to  feed  the  network  proper  and  then 
a  section  requiring  fractional  feed  points  can  be  fed  by  a  system  of  variable 
resistors. 

The  problem  of  properly  designing  a  current  feeding  system  for  the 
network  could  only  be  solved  by  constructing  the  network  of  resistors,  in 
which  no  apparent  difficulty  has  been  found,  and  then  building  up  the  best 
possible  current  feeding  system  and  trying  it  out.  The  requirements  of  this 
current  feeding  equipment  are  now  more  clearly  understood.  Almost  all 
the  time  spent  on  the  simulator  was  spent  in  the  development  of  the  current 
feed  and  it  is  felt  that  a  suitable  system  can  be  devised  as  a  result  of  this 
endeavor.  Section  D  which  follows  shows  some  of  this  work  in  progress, 
and  current  thinking  on  this  problem.  None  of  the  work  can  be  considered 
complete  and  therefore  no  final  design  of  a  cur  rent  feed  system  is  recom¬ 
mended.  in  this  report. 
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SECTION  D 


FURTHER  EXPERIMENTAL  AND  DEVELOPMENT  WORK 


A.  Iteration  Method 


This  method  of  feeding  current  has  been  indicated  in  the  other 
sections  of  this  report.  The  circuit  is  shown  in  fig.  D-l  for  a  two  dimen¬ 
sional  network.  Three  node  points  are  shown  in  fig.  D-l  in  a  network  of 
any  number  (n)  node  points. 

In  Section  A  it  was  shown  that  ii ;  i2  ;  i3  etc.  must  be  equal  to 
A Vq  (l);  AVg  (Z);  A Vg  (3)  etc.  respectively  or 


A  _ii _ h _  __L _  =  _in _ 

"ve(l)  ~  ve  (2)  "  V0  (3)  . . ”ve  (n) 


When  the  scanning  switch  and  amplifier  were  used  the  A  was  the  same  for 
each  feed  point  so  it  was  only  necessary  to  find  the  value  of  A  which  gave 
a  constant  Vg  (l),  Vg  (Z)}  Vg  (3)  etc. 

In  this  method,  however,  the  voltages  are  always  constant  no 
matter  what  the  value  of  A  but  this  current  voltage  ratio  A,  for  every  point, 
must  be  the  same  and  it  is  this  requirement  which  will  be  met  by  only  a 
particular  value  of  A.  The  process  theri  consists  of  finding  this  particular 
A  value  by  adjusting  the  variable  resistors.  As  an  example,  the  method  of 
finding  A  for  a  particular  network  is  shown  in  fig.  D-Z.  The  Iteration 
Process  for  a  particular  network  is  shown  in  fig.  D-Z.  The  network  is 
shown  at  the  top  of  the  figure.  The  total  length  $  has  been  divided  into  ten 
equal  sections  and  the  point  5  is  the  center  feed  point  of  the  network.  The 
purpose  is  to  find  a  value  A  such  that  i  fed  into  every  point  will  be  AVg. 

The  shape  of  the  current  feed  curve  (current  as  a  function  of  £)  will  then 
be  the  same  as  the  shape  of  the  V  g  curve  (Vq  as  a  function  of  jj).  Only 
one  particular  value  of  A  will  satisfy  this  condition.  In  this  case  Step  1 
started  with  a  straight  line  current  distribution.  This  results  in  the  curve 
Fx  (x)  for  Vq  which  is  not  a  straight  line  because  the  end  points  are  at 
ground  or  zero  voltage.  Actually,  it  will  slope  approximately  as  shown  in 
the  Vq  curve  of  Step  1.  Then  Fx  (x)  is  a  nearer  approximation  of  the  re¬ 
quired  current  curve  and  it  is  possible  to  make  i  =  Fx(x)  in  Step  Z  where 
Ax  is  the  current  i  entering  any  network  lattice  unit  divided  by  the  corre¬ 
sponding  voltage  Vg  .  It  would  be  best  to  use  the  current  and  voltage  at 
point  5  because  it  will  be  affected  least  by  subsequent  adjustments,  but 
actually  any  point  may  be  used.  Then  in  Step  Z  a  new  Vg  curve  f2  (x)  is 
obtained. 
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Repeating  this  process  in  Step  3  f3  (x)  is  determined  for  Vj g  which 
is  nearly  one  loop  of  a  sine  curve.  (Theoretically  the  curve  should  be  a 
sine  curve.)  Step  4  makes  the  shape  of  the  curve  nearer  still  to  a  sine 
curve.  The  value  A4,  calculated  after  step  4,  would  be  nearly  equal  to  A3. 
This  means  that  more  steps  will  not  change  the  distribution  of  current 
materially. 

In  general  when  A(n  +  l)  is  nearly  equal  to  A  (n)  the  proper  current 
setting  exists.  This  is  the  value  A  which  is  being  sought. 

It  is  evident  that  had  one  started  with  the  current  distribution 
according  to  fx  (x)  or  f2  (x)  the  iteration  process  would  have  been  shortened. 
This  can  be  done  since  it  is  known  that  the  ends  of  the  final  curve  will 
slope  and  the  peak  will  occur  at  the  center  point  and  there  is  no  reason  for 
starting  with  a  straight  line  current  distribution.  Furthermore,  if  a  point 
of  the  6  network  is  grounded,  which  is  the  case  when  simulating  a  rod, 
black  to  thermal  neutrons,  inserted  at  some  particular  point  in  the  reactor, 
Vq  goes  to  zero  at  that  point.  The  peak  voltage  occurs  at  some  point  be- 
tween  the  perimeter  ground  and  the  rod  and  a  reasonable  distribution  of 
current  can  easily  be  estimated.  This  sort  of  forethought  shortens  the  iter- 
ation  process  for  any  problem.  The  chief  advantage  offered  by  this  method 
is  that  switch  and  amplifier  difficulties  are  eliminated.  This  method  applied 
fo  several  hundred  feed  points  could  prove  tedious,  but  it  is  a  convergent 
process  and  will  eventually  give  an  answer  to  the  problem  which  is  set  up  on 
the  simulator.  In  order  to  test  this  method  of  feeding  current,  an  assembly 
of  13  variable  resistors  was  constructed  and  used  to  feed  the  large  mesh 
half  of  the  two  dimensional  network  (see  fig.  B-9).  The  assembly  is  shown 
on  fig.  D-3.  The  plot  of  voltage  distribution  obtained  with  this  equipment  is 
shown  for  two  radii  on  fig.  D-4  and  fig.  D-5.  The  k  as  determined  from  the 
apparatus  was  1.268  which  compares  well  with  the  analytical  value  1.266. 

If  the  switch  and  amplifier  method  of  feeding  current  eventually 
proves  impractical  this  method  should  be  developed  more  fully.  An  auto- 
matic  electromechanical  system  for  adjusting  the  variable  resistors  can  be 
developed. 


B.  Fractional  Feed  Points 

In  Section  05,  it  was  pointed  out  that  the  current  fed  into  a  network 
point  must  be  proportional  to  the  volume  of  reactor  material  represented 
by  that  point.  In  the  following  discussion  a  definite  volume  of  reactor  is 
considered  a  unit  volume.  The  network  representing  that  area  is  a  unit 
lattice  and  the  current  fed  to  the  node  point  in  the  unit  lattice  is  considered 
a  unit  of  current  feed. 
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Consider  again  the  scanning  switch  described  in  Section  C.  If  each 
segment  represents  a  unit  of  current  feed  (  amplifier  current  x  time  of 
contact)  then  a  lattice  representing  twice  the  unit  volume  can  be  fed  twice 
and  the  current  fulfills  the  requirements  of  the  simulator  design.  How¬ 
ever,  in  the  other  direction  there  is  a  problem.  If  a  volume  of  1/2  is  to  be 
fed  the  segment  must  be  cut  in  half.  Even  if  this  were  practical,  further 
reductions  of  l/4  and  1/8  are  certainly  difficult  to  construct  accurately. 

A  better  method  is  to  feed  all  such  fractional  points  with  another  amplifier 
and  switch.  In  this  case  the  segments  can  be  of  a  set  standard  size  but  the 
amplifier  will  feed  a  current  of  say  l/8  or  l/l6  of  the  current  being  fed  to 
the  main  network.  This  establishes  a  new  system  of  basic  units  which  can 
be  used  in  multiples  of  1,  2,  4  etc.  depending  upon  the  volumes  to  be  repre¬ 
sented. 


Another  way  to  accomplish  this  is  to  use  the  variable  resistor  feed 
system,  described  in  Section  D-A,  to  feed  the  fractional  areas  still  using 
the  amplifier  and  scanning  switch  to  feed  the  rest  of  the  network.  The 
mechanics  of  this  combination  feed  method  have  not  been  worked  out  but 
should  be  tried  in  lieu  of  adding  another  amplifier  and  switch  in  subsequent 
work. 


C.  Re-design  of  Scanning  Switch 

Design  for  an  improved  scanning  switch  is  shown  on  fig.  D-6.  In 
principle  this  switch  inverts  the  action  of  the  switch  in  Section  C.  The 
spring  contacts  are  stationary  and  the  contact  segment  is  moved  over  them. 
Thus  the  contact  time  must  be  more  accurate  because  the  same  segment 
length  is  used  for  each  contact.  In  order  to  test  this  design  a  scale  model 
switch,  having  only  six  contacts  was  built  and  tested.  The  switch  is  shown 
on  fig.  D-7.  The  test  was  set  up  as  shown  in  fig.  D-8. 

The  test  work  showed  that  the  contact  time  was  almost  absolutely 
uniform  for  all  contacts  at  least  as  nearly  as  could  be  measured  with  an 
oscilloscope  screen  picture  of  the  contact  operation.  Figure  D-9  is  the 
photograph  of  the  oscilloscope  screen  picture  of  contact  operation.  To 
make  this  picture  all  six  of  the  contacts  were  wired  together.  Each  time  a 
contact  was  made  the  line  on  the  oscilloscope  screen  was  moved  downward. 
Therefore,  the  short  upper  lines  represent  the  contact  open  time  and  the 
longer  bottom  lines  show  the  contact  performance  of  the  switch.  Where  the 
contact  lower  line  is  clean  and  solid  good  contact  is  indicated.  The  hash  on 
two  of  the  contact  lines  indicates  the  presence  of  dirt  which  caused  a  high 
contact  resistance.  When  this  run  was  started  all  the  contacts  were  clean 
and  none  of  the  contact  lines  on  the  scope  picture  showed  any  hash.  The 
picture  was  taken  after  about  two  hours  of  operation  when  some  dirt  began 
to  accumulate.  The  uniformity  of  contact  was  still  very  good  as  can  be 
seen  from  the  picture.  The  length  of  the  contact  line  represents  time. 
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The  scope  beam  moved  from  left  to  right  and  completed  one  screen  sweep 
in  l/30  sec.  The  contact  time  was  therefore  about  0.003Z  sec. 

The  best  results  (longest  satisfactory  runs)  were  made  with  the 
rotor  very  lightly  oiled  and  continuously  cleaned  with  a  swab.  The  design 
(fig.  D-6)  provided  for  running  the  rotor  in  an  oil  bath  with  a  view  to  re¬ 
ducing  surface  wear  but  the  oil  film  established  a  perfectly  insulated 
surface  on  the  rotor,  and  this  idea  had  to  be  discarded. 

In  addition  to  the  spring  contact  shown  on  the  design  drawing  a 
flat  spring  was  tried,  which  had  no  raised  wire  contact  surface.  This 
contact  did  not  work  as  well  as  the  original  design. 

A  further  idea  was  tried  in  which  a  square  thread  was  cut  on  the 
rotor  so  that  a  groove  was  formed  that  wiped  across  the  stationary  contact 
surfaces  once  in  each  revolution.  This  did  not  offer  any  marked  improve¬ 
ment.  The  contacts  still  became  dirty  after  about  a  six  hour  run. 

No  definite  conclusions  have  been  reached  about  this  switch  design 
up  to  the  present  time  and  it  is  intended  that  several  different  contact 
materials  will  be  tried  in  an  effort  to  find  something  more  satisfactory  than 
steel  on  brass.  These  materials  were  used  only  because  they  were  readily 
available  from  stock.  It  is  felt  that  a  run  of  40  hours  without  cleaning  is 
necessary.  In  general,  the  results  of  test  work  are  encouraging  and  it  is 
felt  that  this  method  of  current  feeding  should  not  yet  be  abandoned,.  Even¬ 
tually  it  may  become  an  established  fact  that  a  sliding  contact  is  not 
suitable  for  this  service.  Other  methods  of  scanning  should  be  investigated, 
especially  an  electronic  multivibrator  type  circuit.  The  question  of  trans¬ 
mitting  a  very  low  voltage  signal  through  such  a  circuit  may  prove  to  be  a 
problem,  since  the  V@  signal  voltage  ranges  from  0.Z  to  Z.O  volt.  The  Vq 
voltage  can,  however,  be  increased  by  changing  the  kappa  factor  of  the 
network  design  and  thus  the  voltage  can  be  raised  to  a  more  suitable  value. 


D.  Re-design  of  Current  Feeding  Amplifier 

In  Section  05  some  of  the  difficulties  encountered  in  the  amplifier 
were  pointed  out.  An  amplifier  which  eliminates  most  of  these  troubles 
was  designed  and  tested  by  F.  W.  Manning  in  the  electronic  shop  at  ORNL». 

The  following  analysis  of  the  problem  and  amplifier  design  was  done  by 
Manning  and  is  included  here  for  completeness. 

It  is  clear  that  this  amplifier  meets  the  specifications  set  forth  by 
the  author  for  the  simulator  amplifier.  While  this  project  was  in  progress 
it  was  decided  to  discontinue  work  on  the  simulator  and  hence  no  performance 
tests  were  made  on  the  network. 


The  reactor  simulator  amplifier  (fig.  D-15)  is  a  direct  coupled  three 
terminal  amplifier  designed  to  take  a  negative  input  signal  of  up  to  three  volts 
and  deliver  a  current  up  to  nine  milliamperes  to  the  output  terminals.  The 
output  current  is  relatively  independent  of  the  load  circuit  impedance.  The 
transconductance  of  the  amplifier  may  be  readily  adjusted  to  any  value  from 
700  to  3000  micromhos.  When  operated  from  a  regulated  a-c  line5  the  trans¬ 
fer  characteristic  of  the  amplifier  remains  essentially  a  constant  negative 
transconductance  at  load  currents  greater  than  50  microamperes. 

The  circuit  for  the  amplifier  consists  of  a  cathode  follower  input 
stage,  a  “bootstrap”  cathode  follower  output  stage,  an  attendant  amplifier, 
and  the  necessary  power  supplies.  The  basic  circuit  is  that  of  a  “bootstrap” 
or,  cathode  follower  used  as  a  constant  current  device,  which  may  be  found 
on  page  94  of  Time  Bases  by  O.  S.  Puckle  (John  Wiley  and  Sons).  This 
circuit  is  shown  in  fig.  D-10. 

In  this  circuit,  the  current  charging  the  capacitor  C  remains  almost 
constant  after  opening  switch  S  until  the  potential  across  the  condenser 
approaches  a  large  fraction  of  the  supply  voltage.  According  to  the  above 
reference,  the  linearity  of  the  charging  current  is  the  same  as  if  condenser  C 
were  being  charged  through  an  ohmic  resistance  from  a  high  potential  supply. 
In  the  case  shown  in  fig.  D-10  with  the  cathode  resistor  made  equal  to  the 
plate  resistance,  the  equivalent  supply  voltage  is: 

ju  +  2  _  1  _ 

- E  +  — Eg 

2  2 

A  modification  of  the  above  circuit,  shown  in  fig.  D-ll  was  utilized 
in  the  Pile  Simulator  Amplifier  fig.  D-15  wherein  a  6AG7,  pentode  connected, 
V6,  preceded  by  a  6SL.7  difference  amplifier,  V8,  is  used  instead  of  the  triode 
above.  A  pentode,  V9,  is  utilized  as  the  cathode  resistor  of  the  “bootstrap” 
which  allows  operation  of  the  “bootstrap”  cathode  follower  with  a  reasonable 
static  plate  current,  thus  operating  in  a  more  linear  portion  of  the  tube 
characteristics,  while  retaining  the  advantage  of  having  a  high  dynamic 
impedance  in  its  cathode  circuit. 

The  input  signal  is  applied  to  Vu,  a  6AG7  connected  as  a  triode  cathode 
follower.  The  output  of  this  stage  is  connected  in  series  with  the  transcon¬ 
ductance  control  R  ,  and  connected  thence  to  the  cathode  of  V8, 

The  bias  control  is  so  adjusted  that  the  static  plate  current  of  Vu 
flows  through  the  bias  control  and  no  potential  drop  occurs  across  Rq. 

The  zero  control  adjusts  the  static  bias  on  V6.  Under  normal  oper¬ 
ating  conditions,  this  control  is  used  to  adjust  the  output  current  to  zero 
when  both  input  and  output  are  shorted. 
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This  control  also  changes  the  effective  gain  of  the  amplifier  Vg  and 
consequently  will  have  a  minor  effect  on  the  over-all  characteristics  of 
the  instrument.  This  is  not  considered  a  serious  handicap  since  the  zero 
point  will  remain  fixed  for  a  given  set  of  tubes  in  the  amplifier.  Zero  drift 
is  minimized  when  the  instrument  is  operated  from  a  regulated  a-c  line. 

Simplification  of  fig.  D-ll  results  in  fig.  D-12  which  lends  itself 
more  readily  to  analysis. 

V6  and  V8  of  fig.  D-ll  are  replaced  by  a  triode  V2  of  fig.  D-12  whose 
plate  resistance  is  identical  to  that  of  Vg  and  whose  amplification  factor  is 
greater  by  the  gain  of  V8.  Vx  of  fig.  D-12  corresponds  to  Vu  of  fig.  D-ll. 
Since  very  little  current  change  occurs  in  V9  little  error  is  encountered  in 
its  omission  in  fig.  D-12.  Since  current  will  be  considered  flowing  into  the 
load,  little  error  is  introduced  by  omitting  Vj0  from  fig.  D-12. 

It  will  be  observed  that  the  change  in  grid  voltage  of  Vx  is: 


degl  =  d6s  "  dipl[RB+RC]  ‘  dip2RC 
and  the  corresponding  change  in  plate  voltage  is: 

depl=_dipl  [RB+RC]  ■  dip2RC 


(1) 


(2). 


By  application  of  the  first  term  of  a  Taylor  series  expansion  for 
plate  current  of  a  triode  vacuum  tube,  the  change  in  plate  current  is. 

di  =  3-2-  de  +  a2E-de  =  g  e  + 

P1  -v  g  x  P  m  g 

de  o  e  rp 

g  P  ^ 


and  for  Vi 


,  der>l 

di  =  g  de  +  — 
pi  “ml  gl 

Pi 


(3) 


From  1,  2  and  3: 


dip  =  gml  [d6s  '  dipl  (RB  +  RC}"  diP2RC]  +  r  F  ^p^W  "  dip2Rc]  (4) 

pi 


P 

Since 


m  r 

P 

dipl  [rp  +  (Kl  +  1)(RB+RC)]  =  Mldes  "  +  RCdip2 


(5) 
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Whence 


Mldes  "  ^1  +  ^RCdli 


pl  rpl  +  (M+  !)(Rb  +  Rc 


Likewise  for  V- 


de  =  -di  (R  +  R  )  -di  _R^, 
g2  p2  A  C  p2  C 


With  a  load  resistance  R  across  the  output  terminals 

i_j 


de  =  -di  (R  +  R  +  R  )  -di  R 
p2  p2  L  A  C  pl  C 


di  =  g  de  +  — 
p2  mZ  gZ  r 


diP2  =  8m2[-dip2  (Ra  +  Rc)  -diplRc]  +  ±  [-dip2(RL  +  Ra  +  Rc)  -diplRc]  (10) 
V  [rp2  +  +  1,(RA  +  Rc)  +  Rl]  =  -(^2  +  l)RCdipl  (11) 


Substituting  (6)  in  (ll) 


V  [[rp2  +<»*2  +  1)(RA  +  RC)  +  RJ  ' 


(mx  +  i)(m2  +  i)r; 


rPi+(Mi  +  lKRB  +  Rc 


rcMi<p2  +  i) 


r  +(Ml  +  1)(R  +  R  )  s 


Since  the  transconductance  of  a  network  is  defined  as  the  ratio  of 
the  change  in  output  current  to  the  change  in  input  voltage,  from  12  above 
it  is  found  that: 

%  M  |  (  M  2  +  1  )Rq  (13> 

A  =  =  ’  t  VU1  +  1,(RB+RC>]  [  rp2  +(fl2  +  1)(RA+  RC)  +  RJ  ‘[^l  +  ‘Mm^IJR; 
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By  substituting  the  following  values  for  the  quantities  as  found  in 
reactor  simulator  amplifier: 

r  ,  =  2500 
P1 


Rb  =  150 

Rc  =  2200 

r  =  130,000 

p2 

u  =  Gg  r  =  30  x  11000  x  10"6  x  130,000  = 

p2  sm  p 

Ra  =  Ra 
A  A 


42800 


Equation  (13)  then  reduces  to: 


A  =  - 


1.172  R  +  282 

A 


Since  a  1000  ohm  10  turn  helipot  was  used  for  R  ,  with  the  dial 
calibrated  from  0  to  1000,  but  so  oriented  that  the  resistance  in  the  circuit 
was  equal  to  1000  minus  the  dial  reading,  the  above  equation  reduces  to  the 
following  in  terms  of  dial  reading. 


A  =  - 


1454  -  1.172D 


A  plot  of  this  equation  and  of  experimental  data  taken  with  a  one 
volt  signal  applied  to  the  input  terminals  and  a  100,000  ohm  resistor  across 
the  output  terminals  is  shown  in  fig.  D-13. 

It  will  be  noted  that  close  agreement  exists  between  the  two  curves, 
substantiating  the  above  circuit  analysis. 

Typical  input  voltage  -  output  current  curves  are  shown  in  fig.  D-14 
for  two  settings  of  the  trans conductance  control.  It  may  be  seen  that  these 
curves  are  essentially  straight  lines.  The  curve  for  a  dial  setting  of  975  with 
a  100,000  ohm  load  illustrates  the  saturation  of  the  amplifier  at  an  output 
voltage  of  420  volts  corresponding  to  4.2  milliamperes  through  the  load 
res  istance. 

The  action  of  the  amplifier  operating  into  a  capacitive  load  may  be 
analyzed  by  assuming  a  capacitor  C  across  the  output  terminals.  The  dif¬ 
ferential  equation  for  the  current  flowing  in  the  output  terminals  may  be 
shown  to  be: 
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(/Li1  +  l)(jU2+l)  R 


M  1(^2+1)Ri 


rPi+  c  s 


*p2  rp2+^M2+1^RC+RA^ 


which  is  of  the  form 
IR+§=E 


with  solution 


1  =fe  "  Rc  (18) 

•  From  this  it  may  be  seen  that  the  capacitor  C  is  apparently  being 
charged  from  a  voltage  source 

n  fn  I  1  ^  R  (^) 

~  _  ^1^2  }  C _  _  21  x  4Z800  x  2200  _ 

'  "  es  r  x  +  (/x1+l)(R  +Rb)  "  "  es  2500  +  22(2200  +  150)  “  “  36600  es 


The  apparent  series  resistance  R  will  be  found  to  be 


(M1  +  1)(m2+1)  Rc2 


*-pa  +  (yi)(w-  r;1,(,1;.)(Rc,HB) 


1.03  x  107  +  42800  R 


Thus  the  apparent  series  resistance  will  vary  from  about  10  megohms 
to  53  megohms. 


Experimental  results  showed  approximately  10  microamperes  change 
in  output  current  by  shorting  out  a  100,000  ohm  load  resistor,  while  the  trans1 
conductance  was  adjusted  to  ZOOO  micromhos.  This  implies  an  apparent 
voltage  of  approximately  -40,000e  and  apparent  series  resistance  of  about 
Z0  megohms.  This  shows  a  reasonable  correlation  between  calculated  and  ex1 
perimental  results. 

•The  apparent  time  constant  of  the  circuit  when  used  with  a  load 
capacitance  of  10  microfarads  would  be  of  the  order  of  100  to  530  seconds. 
For  periods  small  with  respect  to  the  time  constant  one  would  expect  to 
charge  the  condenser  with  a  reasonably  constant  current. 


It  was  found  necessary  to  operate  the  amplifier  from  a  regulated  a-c 
line  to  prevent  excessive  drift  of  the  output  current.  This  occurred  in  spite 
of  the  fact  that  an  electronically  regulated  supply  with  a  stabilization  factor 
of  1:200  was  used  to  furnish  the  positive  and  negative  105  volts  to  the  am¬ 
plifier  stage,  V8  of  fig.  D-ll.  It  is  feit  that  considerable  benefit  would  be 
obtained  by  balancing  this  stage  by  the  addition  of  a  one  megohm  plate  load 
to  the  input  section  and  by  adding  a  15  megohm  resistor  from  its  plate  to 
negative  105  volts.  This  would  tend  to  prevent  changes  in  supply  voltage 
from  affecting  the  bias  of  the  “bootstrap”  stage.  This  would  produce  an 
attendant  loss  in  gain  and  consequently  reduces  the  value  ofp2  by  approxi¬ 
mately  l/2,  thereby  modifying  all  of  the  above  calculated  and  experimental 
results.  The  input  impedance  of  the  amplifier  is  approximately  equal  to 
the  input  grid  resistor  since  the  signal  is  almost  completely  degenerated  in 
the  cathode  circuit  of  Vn.  Should  a  higher  input  impedance  be  desired,  it 
is  suggested  that  Vu  be  preceded  by  a  high  input  impedance  cathode  follower 


^Master 
Voftaye  Coniro/ 
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cy 


p.pa/4/p  Fiq-D't 


X  U 


Network  L/ 


Barren] 


BBS 


/  2.  3  4  *  a  7  a 


Const. 


JL  - 

Vo  -f,  ex) 


4  5  6  7  8  9  to  O  /  Z  3  4  5  6  7  3  9  /0 

— *■  X 


L  :Aifi  (X) 


Vo'-TzM 


C=  Aih  (X) 


V*  =  *  (X) 


C  *  As  is  60 


When  A  »  A„  Then  The  Critical  A  Has  Been  Found 

h+f  ft 


Iteration  Process 
p.qd.  4-/46  Fig.D-Z 


PLOT  OF  <pf  AND  <£th  FOR  SIMULATOR  ADJUSTED  BY  ITERATION 
METHOD.  DISTRIBUTION  ALONG  RADIUS  OF  SIMULATOR  THROUGH 
,00  _  POINTS  01;  022;  044;  066. 
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S110A 


PPD  -  A  -  1399 


METHOD.  DISTRIBUTION  ALONG  CHORD  OF  RE  ACTOR  -  S I M  UL  ATOR 
IOO  -  POINTS  01;  020;  040;  060;  080;  0100. 

* -  REACTOR  REGION - - dm  REFLECTOR  REGION - J 
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sitoa 


FIG.  D-  5  PPD-A-  1400 


TRANSCONDUCTANCE  -  MICROMHOS . 
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Fig.  D-15 


J-ft  T-4  Sioncar  Transformer  P40J6,  325-Q-J25V 
@>40ma.,  5. 0  V.  <®  2  0 Amp-,  6.3  V.  @  20  Amp 

T-2,  7-  3/  Thondarson  Transformer  -  T2./F/0 
63  V(8>  3  Amp 

T~S  A/ tied  '  /.-/  I  so/a //an  Transformer  if '137 
HSH-USV.  p  30-  40 Y A.  5^' 


APPENDIX  I 


Introduction 


In  the  body  of  this  report  it  is  shown  that  practical  construction 
considerations  necessitate  the  use  of  a  lumped  network  for  the  electrical 
analogy  of  a  neutron  reactor.  It  was  pointed  out  by  J.  W.  Simpson  that  the 
differential  equation  of  an  electrical  transmission  line,  with  distributed 
sources  and  sinks,  was  exactly  similar  to  the  differential  equation  of  a 
slab  reactor  in  elementary  diffusion  theory.  Having  drawn  this  analogy  one 
then  can  build  a  "tt”  section  of  resistors  which  is  exactly  equivalent,  at  its 
terminals,  to  the  transmission  line.  A  lumped  network  of  these  “  7T  sec¬ 
tions  would  then  be  exactly  equivalent  to  a  distributed  slab  reactor  in  ele¬ 
mentary  diffusion  theory.  The  formulas  for  the  equivalent  circuit  of  a 
transmission  line  are  well  known  and  can  be  found  in  any  text  book  on  the 
subject.  Here  the  necessary  formulas  for  such  a  resistor  network  have 
been  derived  and  the  analogy  extended  to  include  the  two  group  theory  of 
neutron  diffusion.  The  formulas  for  a  three  dimensional  reactor  are  also 
derived. 

The  formulas  for  the  relation  between  the  various  electrical  network 
constants  involve  the  value  k  of  the  reactor  to  be  simulated.  This  means 
that  the  analogy  would  only  be  exact  for  one  k  value.  For  any  value  of  k 
other  than  this  particular  value,  the  relations  are  again  approximate.  The 
details  of  this  analogy  are  included  here  for  completeness  and  because 
these  more  exact  formulas  may  prove  useful  in  the  future.  No  use  was 
made  of  these  formulas  in  the  experimental  work  to  date. 


APPENDIX  I 


Transmission  Line  Analogy 


The  results  of  the  finite  difference  calculation  can  be  reached  in  a 
more  exact  form  by  making  the  following  analogy  between  a  transmission 
line  and  a  reactor.  It  is  possible  to  build  lumped  network  exactly  equiv¬ 
alent  to  a  distributed  transmission  line  at  the  end  points  of  the  line.  By 
making  the  direct  analogy,  then  one  can  build  a  network  exactly  equivalent 
to  a  distributed  reactor. 


A.  One-Group  Theory 


tor. 


Consider  a  slab  pile  of  any  length  l  surrounded  by  a  perfect  reflec- 


s> 


g 


Neutron  Flow 


Moderator" 


4> 


r 


i 


(0)  is  the  flux  at  a  distance  1  from  (4>)r. 


Fissionable  material  is  distributed  homogeneously  along  the  reactor  mod¬ 
erator. 


Reactor  material  constants  are  Za,  3Z^ 

K  =  number  of  thermal  neutrons  produced  per  thermal  neutron  absorbed  in 
an  infinite  reactor. 

The  equivalent  electrical  transmission  line  is* 
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In  which  volt  g  is  the  current  generated  per  volt  every  dl  distance  along 
the  line  and  is  analogous  to  in  the  reactor. 

y'=  2  a  mohs  per  unit  length. 

r  =  ohms  per  unit  length. 

Then 

I  =  Current  in  amps,  is  analogous  to  neutrons  flowing  at  any  point. 

E  =  Voltage  anywhere  along  the  line. 

Er  =  At  receiving  end  =  (nv)r  =  0  r 
Eg  =  At  sending  end  =  (nv)g  -  <J>g 

I  =  Length  (cm)  and  will  be  measured  from  the  receiving  end. 


dE 

df 


rl 


d2E  rdl 
d£2  “  df 


Let  the  net  conductance  y  =  y'-  g  (since  g  is  a  negative  conductance). 
Then: 


d2I  _  dE 
df 2  ~  7  d i 


d2E 

d£2  " 


ryE 


d21 

dF  =  ryI 


(1) 


Letting  D  = 


di 


the  solution  of  these  equations  is 


(D2  -  ry)  E  =  O 


(D2  -  ry)  1=0 


(2) 


E  =  A  exp  ft  \fry  +  B  exp  -  ft  Vry 
I  =  C  expi  Vry  +  D  exp  -  ft  Vry 
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With  receiving  end  short  circuited.  Er  =  O 

Eg  =  I rJy  sinh  {  \/ry  (7) 

Ig  =  Ir  cosh  P  \fry  (®) 

A  lumped  resistance  network  which  is  equivalent  to  the  line  of  dis¬ 
tributed  constants  must  be  found.  This  can  be  obtained  in  the  following  way: 


Let  Rf  be  the  resistance  to  ground  from  (l)  and  (Z)  in  (a) 
Let  Rg  be  the  resistance  to  ground  from  (l)  in  (b) 

Then 


Rf  =  Rx  +  R3 

O  -  P  .  R1R3  Rl  2  +  R1R3  +  R1R3 

Kg  Kl  Rj  +  R3  Rj  -  R3 

Eliminate  Rj 


Solve  for  R3  from  (7)  and  (8) 

R  =  _^g  tanh  A  \/ry 

5  T  y 
xg 

Likewise  from  (5)  and  (6) 
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Insert 

y  =  y'  -  g 

The  general  relation  is  obtained 

—  -  RG  =  2  tanh  \/~r  (y'~  g)  ,  sinh  l  \J  r  (y'  -  g)  (10) 

R2  * 

in  which  G  and  g  are  any  negative  conductances,  the  values  of  which  can  be 
specified  as  desired. 


The  equivalent  7 r  section  will  be: 


from  which  it  is  evident  that  the  same  relations  will  hold  for  the  resist 
ances. 


This  method  can  be  used  for  three  dimensions  with  equal  results 
by  simply  considering  that  the  neutrons  flow  only  in  the  directions  x,  y, 
and  z.  Each  direction  is  then  a  separate  transmission  line  having  the  re 
sistance  relations  of  formula  (10). 


B.  Two-Group  Theory 

This  method  is  now  extended  into  two-group  theory  by  using  the 
results  of  the  preceding  calculations. 


For  simplicity  a  slab  pile  is  again  considered  and  an  equivalent 
electrical  transmission  line  is  drawn  for  each  neutron  energy  group. 


High  Energy  Group 


Low  Energy  Group 
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1 


The  high  energy  group  line  is  put  on  top  of  the  low  energy  group 
line  in  such  a  way  that  E^2af  no  longer  feeds  to  ground  but  instead  feeds 
into  the  low  energy  line  in  place  of  gx  as  follows: 


w 


3Z^f  per  length 

-AAAA/v — 

32  per  length 

— wwv — 


We  are  interested  in  finding  such  values  for  y 0  and  y^  that  with 
the  lines  combined,  each  will  have  the  same  current  flow  as  it  did  with  the 
lines  separated. 


For  the  separate  lines  if  g  =  h2ath  following  equations  for 
current  are  obtained. 


~ af  =  E0  2af"k2athE@  (11) 

and 

-di-=2athE0  “  2af  E0  (12) 

For  the  combined  lines 

-Tf  =  (E0  ‘  V  y*  -  sEe 

■Tf  =  E*  ye  -(s+Yt)  Eg  (13) 

and 

^dl"  =  E@  ye  -  (E0  -  Eg)  y<p 

■^af  =  <ye  -y</>)  Ee "  y0 E0  (14) 

So  that  (ll)  and  (12)  will  equal  (13)  and  (14)  respectively,  the  following 
conditions  must  exist: 

V(j>  =  2af 
y@  ~  ^ath  af 


(15) 


g'  =  2  af 


(IS) 


These  relations  will  be  used  later. 

A  T  network  of  resistors  can  be  treated  in  a  similar  manner.  The 
T  sections  are  set  up  as  follows: 


Put  the  0  network  on  top  of  the  Q  network 


R q  ,  and  R-y  are  the  same  as  before  and  Rg  is  an  adjusted  value  to  be 

determined.  The  currents  V ^  l'Q  are  those  resulting  from  the  rearrange 
ment  of  the  networks. 
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Since  the  currents  in  the  resistors  R0  and  Rq  are  to  be  the  same 
in  the  combined  networks  as  in  the  separate  networks,  l'a  is  less  by 


ifl  = 


r  -  iEs  =  Ei 
e  -  Rg  -  r-  -  r7 


R  g  Ry  Rg 

G,  for  the  networks  separately  can  be  defined 

_  _k_  _  k_  k_ 

1  =  Rg=Rg"R7 

This  conductance  must  be  reduced  by  *pr  mhos  so  that  the  currents  in  Rg 
and  R0  will  not  change,  and,  as  a  result/ the  conductance  for  the  combined 
networks  is 

”  Rg  ‘  R7  "  Ry  (17) 

Using  equation  (10X)  to  relate  the  distributed  values  to  the  lumped  values 
the  following  are  obtained  for  separate  networks  and  lines: 

For  0  Network 

^  -  R0  G  =  2  tanh  £  ^32^2^  -  32tf  g  sinh*  ^32tf2af  -  32^T 
For  8  Network 

Rg  '  R0  Gi  =  2  tanh-A\/32tth2ath  -  32  tth  g  .  smh  J?  \/32tth2ath  "  32tth 

Using  relations  11  through  17  these  can  be  converted  into  similar  relations 
for  combined  networks  and  lines. 


1=2  tanhi-V|-  3kZtf2ath  .  Binltf  $  -  3k2tf2ath 

(18) 


For  0  Network  portion 


1 


For  0  Network  Portion 

=  2  tanh  -  3Ztth2af  .  sinhjj/  ^  -  32tth2af  (19) 

8 

For  small  angles  tanh  7  =  sinh  7=7  and  using  this,  reduction  of  (18)  will 
give  formula  (10)  of  Section  A. 

M 

R7  "  T 

and  for  (19)  formula  (14)  of  section  A  is  obtained. 

Rg  Re  _  A2 
Rg  “  R7  LA 

To  apply  formulas  (18)  and  (19)  to  a  three  dimensional  network  it 
is  necessary  to  consider  that  the  neutrons  move  in  x,  y  and  z  directions  of 
the  cartesian  coordinates.  The  analogy  now  calls  for  three  transmission 
lines  which  meet  at  a  point  x,  y,  z. 


If  each  line  is  then  treated  alone  the  formulas  18  and  19  will  apply 
for  each  direction.  If  a,  b,  and  c  are  made  the  lengths  in  the  x,  y  and  z 
direction  respectively  we  will  get: 


